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ABSTRACT 

We derive an accurate mass distribution of the rich galaxy cluster C10024+1654 (z = 0.395) based on deep Sub- 
aru BRcz' imaging and our recent comprehensive strong lensing analysis of HST/ACS/NIC3 observations. We 
obtain the weak lensing distortion and magnification of undiluted samples of red and blue background galaxies 
by carefully combining all color and positional information. Unlike previous work, the weak and strong lensing 
are in excellent agreement where the data overlap. The joint mass profile continuously steepens out to the virial 
radius with only a minor contribution ^ 10% in the mass from known subcluster at a projected distance of 
~ 700 kpc/;"'. The cluster light profile closely resembles the mass profile, and our model-independent M/L« 
profile shows an overall flat behavior with a mean of (M/Lr) ~ 21>Qh(M /Lr)q, but exhibits a mild declining 
trend with increasing radius at cluster outskirts, r ^ 0.6rvir. The projected mass distribution for the entire clus- 
ter is well fitted with a single Navarro-Frenk- White model with a virial mass, Mvh- = (L2 ± 0.2) x 1O'^M0 /i"', 
and a concentration, Cvii = 9.2^\ 2- This model fit is fully consistent with the depletion of the red background 
counts, providing independent confirmation. Careful examination and interpretation of X-ray and dynamical 
data, based on recent high-resolution cluster collision simulations, strongly suggest that this cluster system is 
in a post collision state, which we show is consistent with our well-defined mass profile for a major merger oc- 
curring along the line of sight, viewed approximately 2 - 3 Gyr after impact when the gravitational potential has 
had time to relax in the center, before the gas has recovered and before the outskirts are fully virialized. Finally, 
our full lensing analysis provides a model-independent constraint of M2d(< '"vir) = (1.4 ±0.3) x 1O'^M0 /?"' for 
the projected mass of the whole system, including any currently unbound material beyond the virial radius, 
which can constrain the sum of the two pre-merger cluster masses when designing simulations to explore this 
system. 

Subject headings: cosmology: observations — galaxies: clusters: individual (C10024H- 1654) — gravitational 
lensing 



1. EVITRODUCTION 

C10024H-1654 (z = 0.395) is the most distant cluster of 
galaxies discovered by Zwicky (1959), and is the focus of 
some of the most thorough studies of cluster properties, in- 
cluding the internal dynamical (Czoske et al. 2002; Diaferio 
et al. 2005), X-ray emission (Soucail et al. 2000; Ota et al. 
2004; Zhang et al. 2005), and both weak (Kneib et al. 2003; 
Hoekstra 2007; Jee et al. 2007) and strong (Colley et al. 1996; 
Tyson et al. 1998; Broadhurst et al. 2000; Comerford et al. 
2006; Zitrin et al. 2009b) lensing work. 

Despite the round and concentrated appearance of this clus- 
ter, several independent lines of evidence point to recent 
merging of a substantial substructure. The internal dynamics 
of about 300 spectroscopically-measured galaxies has been 
modeled by a high-speed, line-of-sight collision of two sys- 
tems with a mass ratio of the order of 2:1, leading to a 
compressed distribution of velocities along the line of sight 
(Czoske et al. 2001, 2002). A direct line-of-sight merger is 
also used to account for the "ring" of dark matter claimed 
by Jee et al. (2007), based on the central mass distribu- 
tion derived from deep Hubble Space Telescope (HST) Ad- 
vanced Camera for Surveys (ACS) images. On a larger scale, 
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the mass distribution derived from a mosaic of HSTAVFPC2 
pointings (Kneib et al. 2003) reveals an additional substantial 
subcluster in the mass at a projected radius of 3' coincident 
with a noticeable concentration of galaxies. This substructure 
is however not along the line of sight, but is associated with 
the main cluster component in redshift space (Czoske et al. 
2002), and lies R ^ 700kpc/!"' northwest (NW) in projection 
from the center of the main cluster. 

No clear evidence of excess X-ray emission is found at the 
location of the NW galaxy clump, but interaction may be im- 
plied instead by the anomalously low level of X-ray emis- 
sion relative to the standard X-ray luminosity-mass relation. 
The measured gas temperature is also unusually low, only 
Tx ~ 4.5 keV (Ota et al. 2004), over the full range of radius 
(R ^ 300kpc/z"') probed by deep Chandra data, whereas re- 
cent weak and strong lensing observations indicate that the 
cluster is a high-mass system with a total projected mass 
of M2D ^ IO'^Mq (Hoekstra 2007; Jee et al. 2007; Mahdavi 
et al. 2008; Zitrin et al. 2009b). The careful hydrodynami- 
cal simulations of Ricker & Sarazin (2001) predict that in a 
period of 1-3 Gyr (or a timescale of the sound crossing time, 
Tse) after a substantial merger the hot gas associated with the 
whole system is extensively distributed so that the emissiv- 
ity is actually markedly reduced by virtue of the lowered gas 
density, once the shock associated with the collision has dis- 
sipated. Furthermore, substantial proportion of the gas may 
escape from the system together with high velocity galaxies, 
so that the velocity dispersion of the remainder is significantly 
reduced. These lack of any observed hot shocked gas compo- 
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nent implies very clearly that no merger has happened very re- 
cently, within a couple of Gyr, unlike the bullet cluster (Clowe 
et al. 2006) and other clusters caught in the first colUsional en- 
counter (Okabe & Umetsu 2008). 

C10024-I-1654 displays one of the finest examples of grav- 
itational lensing forming a symmetric 5-image system of a 
well-resolved galaxy, which was first noted by Koo (1988) 
and later resolved into a close triplet of arcs by Kassiola 
et al. (1992), with two additional images found by Small 
et al. (1996) and by CoUey et al. (1996) using HST WFPCl 
and WFPC2 data, respectively. These arcs have been used 
by Colley et al. (1996) to construct an image of the source, 
whose redshift z = 1 .675 (Broadhurst et al. 2000) permits 
an accurate and model-independent enclosed mass for the 
central R < lOOkpcr^ area of M(< R) = (1.11 ± 0.03) x 
lO^'^M©/! with a central mass-to-light ratio of (M/L)b = 
(32O±3O)/z(Mo/L0)b (Broadhurst et al. 2000). 

More recently the central mass profile has been constrained 
by lensing with the identification of many new multiply lensed 
images (Zitrin et al. 2009b) in very deep multi-color imaging 
with HST/ACS as part of the ACS/GTO program (Ford et al. 
2003). Here a joint fit was made to 33 lensed images and their 
photometric redshifts with a relatively simple 6-parameter de- 
scription of the deflection field. This modeling method has 
been recently applied to two unique, X-ray luminous high- 
z clusters, MACS J1149.5H-2223 and MACS J0717.5-I-3745, 
uncovering many sets of multiply-lensed images (Zitrin & 
Broadhurst 2009; Zitrin et al. 2009a). The high resolution 
and accuracy of colors allow a secure identification of counter 
images by delensing the pixels of candidate lensed images to 
form a source which is then relensed to predict the detailed 
appearance of counter images, a technique developed for sim- 
ilar high quality ACS/GTO data of A 1689 (Broadhurst et al. 
2005a). This model has been used to estimate the magnifica- 
tion of high-redshift candidate galaxies with photometric red- 
shifts Zphot ~ 6-7 identified in combining this data with deep 
near-infrared images (Zheng et al. 2009). Here we add to this 
new strong-lensing information with new weak-lensing data 
from deep, multi-color imaging with the Subaru telescope to 
examine the mass distribution in detail over the full profile of 
the cluster. These high-quality Subaru images span the full 
optical range in the B, R^, and z' passbands, allowing us to de- 
fine a background sample of red and blue galaxies free from 
cluster member and foreground galaxies. We have learned 
from our earUer work that without adequate color informa- 
tion, the weak-lensing signal can be heavily diluted particu- 
larly towards the cluster center by the presence of unlensed 
cluster members (Broadhurst et al. 2005b; Medezinski et al. 
2007; Umetsu & Broadhurst 2008), so that weak lensing un- 
derpredicts the Einstein radius derived from strong-lensing 
studies and the gradient of the inner mass profiles based on 
weak lensing is also underestimated. Unfortunately, many ex- 
amples of this problem are present in the literature, and here 
we carefully explore the weak-lensing signal in color-color 
space and by comparison with the deep photometric redshift 
survey in the COSMOS field (Ilbert et al. 2009). 

A major motivation for pursuing improved lensing mea- 
surements is the increased precision of model predictions for 
the mass density profiles of cluster-size massive dark-matter 
halos based on A^-body simulations in the standard A cold 
dark-matter (hereafter ACDM) model (Hennawi et al. 2007; 
Neto et al. 2007; Duffy et al. 2008). Clusters of galaxies 
provide a definitive test of the standard structure-formation 
model because their mass density profiles, unlike galaxies, are 



not expected to be significantly affected by cooling of baryons 
(e.g., Blumenthal et al. 1986; Broadhurst & Barkana 2008). 
This is because the high temperature and low density of the 
intra-cluster gas (hereafter ICG) prevents efficient cooUng and 
hence the majority of baryons simply trace the gravitational 
potential of the dominant dark matter. Massive clusters are 
of particular interest in the context of this model, because 
they are predicted to have a distinctively shallow mass pro- 
file (or low concentration) described by the form proposed by 
Navarro et al. (1997) and this question has been the focus of 
our preceding work (Broadhurst et al. 2005b; Umetsu et al. 
2007; Broadhurst et al. 2008; Umetsu & Broadhurst 2008; 
Umetsu et al. 2009). 

The paper is organized as follows. We briefly summarize 
in §2 the basis of cluster weak gravitational lensing. In §3 
we describe the observations, the photometry procedure, the 
sample selection, and the weak-lensing shape analysis. In §4 
we present our weak lensing methods, and derive the clus- 
ter lensing distortion and convergence profiles from Subaru 
weak lensing data. In §5 we examine in detail the cluster 
mass and light profiles based on the joint weak and strong 
lensing analysis. In §6 we compare our results with pre- 
vious studies of C10024H-1654 to examine the long-standing 
puzzle on large mass discrepancies between lensing and X- 
ray/dynamical methods, and investigate the implications of 
observed discrepancies and anomalies. In §7 we explore and 
discuss a possible interpretation of the observed X-ray fea- 
tures and mass discrepancies Finally, a summary is given in 
§8. 

Throughout this paper, we use the AB magnitude system, 
and adopt a concordance ACDM cosmology with fi^ = 0.3, 
Oa = 0.7, and h = //o/(100kms-' Mpc"') = 0.7. In this cos- 
mology, 1' corresponds to 224 kpc (and 1" to 3.73 kpc /z"') 
at the cluster redshift. All quoted errors are 68.3% con- 
fidence limits unless otherwise stated. The reference sky 
position is the center of the central bright elliptical galaxy 
(the galaxy 374 in the spectroscopic catalog of Czoske et al. 
(2002)): R.A. = 00:26:35.69, Decl. = h-17:09:43.12 (J2000.0). 
We refer to this position as the optical cluster center, here- 
after. The cluster center of mass for our radial profile analysis 
is chosen to be the dark-matter center at AR.A. = -2.32", 
ADecl. = -1.44" of Zitrin et al. (2009b). 

2. BASIS OF CLUSTER WEAK LENSING 

Weak gravitational lensing is responsible for the weak 
shape distortion and magnification of the images of back- 
ground sources due to the gravitational field of intervening 
foreground clusters of galaxies and large scale structures in 
the universe (e.g., Umetsu et al. 1999; Bartelmann & Schnei- 
der 2001). The deformation of the image can be described by 
the 2 X 2 Jacobian matrix Aap (a, ^ = 1,2) of the lens map- 
ping.^ The Jacobian Aa^ is real and symmetric, so that it can 
be decomposed as 

■A.a0 = {'^-n)^al3-'^a0, (1) 

where daji is Kronecker's delta, T^p is the trace-free, sym- 
metric shear matrix with 7^ being the components of spin- 

^ Throughout the paper we assume in our weak lensing analysis that the 
angular size of background galaxy images is sufficiently small compared to 
the scale over which the underlying lensing fields vary, so that the higher- 
order weak lensing effects, such as flexion, can be safely neglected; see, e.g., 
Goldberg & Bacon (2005); Okura et al. (2007, 2008). 
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2 complex gravitational shear 7 := 71 + ij2, describing the 
anisotropic shape distortion, and k is the lensing convergence 
responsible for the trace-part of the Jacobian matrix, describ- 
ing the isotropic area distortion. In the weak lensing limit 
where k, I7I <C 1, Tap induces a quadrupole anisotropy of the 
background image, which can be observed from ellipticities 
of background galaxy images. The flux magnification due to 
gravitational lensing is given by the inverse Jacobian determi- 
nant, 

l^=-m = - — r-n> (3) 



det^ 



I7|- 



where we assume subcritical lensing, i.e., dttAiff) > 0. 

The lensing convergence is expressed as a line-of-sight pro- 
jection of the matter density contrast 5m = {p„-p)/p out to 
the source plane {s) weighted by certain combination of co- 
moving angular diameter distances r (e.g., Jain et al. 2000), 



2c2 

Q{x,Xs) = 



7 

Jo 



dxQix,Xs)^ = Jd^„J:-^, (4) 



Kx)Kx.-x) 

riXs) 



(5) 



where a is the cosmic scale factor, x is the co-moving 
distance, S„ is the surface mass density of matter, E„ = 
/ dxa{pm—p), with respect to the cosmic mean density p, and 
Edit is the critical surface mass density for gravitational lens- 
ing, 

Scrit = -j— ^ „ n 

A-kG DdDds 

with Dj, Dd, and Djs being the (proper) angular diameter dis- 
tances from the observer to the source, from the observer to 
the deflecting lens, and from the lens to the source, respec- 
tively. For a fixed background cosmology and a lens redshift 
Zd. 5]crit is a function of background source redshift Zs- For 
a given mass distribution Em(^), the lensing signal is propor- 
tional to the angular diameter distance ratio. 



l^izs) = max 



DdsiZs) 
DsiZs) 



(7) 



where p{zs) is zero for unlensed objects with Zs < Zd- 

In the present weak lensing study we aim to reconstruct 
the dimensionless surface mass density k(9) from weak lens- 
ing distortion data. For a two-dimensional mass reconstruc- 
tion, we utilize the relation between the gradients of k and 7 
(Kaiser 1995; Crittenden et al. 2002), 



Ak(6>) = d'^d^T^^iO) = 2©*7(6») 



(8) 



where V is the complex differential operator V = {dl- d^) /2+ 
id\d2- The Green's function for the two-dimensional Pois- 
son equation is A"' (9,6') = ln\9-9'\/ (Itt), so that equation 
(8) can be solved to yield the following non-local relation be- 
tween K and 7 (Kaiser & Squires 1993): 

k(6») = ^jd^9'D*ie- e>')7(6'') (9) 

where D{9) is the complex kernel defined as D{9) = {9\ -9\- 
2i9\92)l\9\'^ . In general, the observable quantity is not the 
gravitational shear 7 but the complex reduced shear. 



g = 



7 

\-K 



(10) 



in the subcritical regime where det^ > (or 1 / g* in the neg- 
ative parity region with det^ < 0). We see that the reduced 
shear g is invariant under the following global transformation: 



k(0)^ Ak(6')+1-A, 7(6»)^A7(0) 



(11) 



with an arbitrary scalar constant A 7^ (Schneider & Seitz 
1995). This transformation is equivalent to scaling the Ja- 
cobian matrix A(9) with A, A{0) \A{9). This mass-sheet 
degeneracy can be unambiguously broken by measuring the 
magnification effects, because the magnification fi transforms 
under the invariance transformation (1 1) as 



p{9) ^ X^p{9). 



3. SUBARU DATA AND ANALYSIS 
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In this section we present a technical description of our 
weak lensing analysis of C10024-H1654 based on deep Subaru 
BR^z' images. The reader only interested in the main result 
may skip directly to §5. 

3.1. Subaru Data and Photometry 

For our weak-lensing analysis of C10024H-1654 we retrieved 
from the Subaru archive, SMOKA, ^ imaging data in B, R^, 
and z' taken with the wide-field camera Suprime-Cam (34' x 
27'; Miyazaki et al. 2002) at the prime-focus of the 8.3m Sub- 
aru telescope. The cluster was observed in the course of the 
PISCES program (Kodama et al. 2005; Umetsu et al. 2005; 
Tanaka et al. 2005). The FWHM in the co-added mosaic im- 
age is 1.27" in B, 0.80" in R^, and 0.82" in z' with 0.202" 
pixer\ covering a field of ~ 34' x 26'. The observation de- 
tails of C10024H-1654 are listed in Table 2. We use the R^-hand 
data for our weak lensing shape measurements (described in 
§3.2) for which the instrumental response, sky background 
and seeing conspire to provide the best-quality images. The 
standard pipeline reduction software for Suprime-Cam (SD- 
FRED, see Yagi et al. 2002; Ouchi et al. 2004) is used for flat- 
fielding, instrumental distortion correction, differential refrac- 
tion, sky subtraction, and stacking. Photometric catalogs are 
constructed from stacked and matched images using SExtrac- 
tor (Bertin & Arnouts 1996). A composite BR^z' color image 
of the central 8' x 8' region of the cluster is shown in Figure 1 . 
Since our lensing work relies much on the colors of galaxies, 
special care has to be paid to the measurement of colors from 
BR^z' images with different seeing conditions. For an accurate 
measurement of colors the Colorpro (Coe et al. 2006) routine 
is used; this allows us to use SExtractor's AUTO magnitudes 
for total magnitudes of galaxies and its isophotal magnitudes 
for estimation of colors, and applies a further correction for 
the different seeing between different bands. From the colors 
and magnitudes of galaxies we can safely select background 
galaxies (see §3.3) for the weak-lensing analysis. 

Astrometric correction is done with the SCAMP tool 
(Bertin 2006) using reference objects in the NOMAD cat- 
alog (Zacharias et al. 2004). Photometric zero-points were 
calculated from associated standard star observations taken 
on the same night. Since for B and Rc only one standard 
field was taken, and for z'-band only a few spectrophotomet- 
ric standards are available, we further calculated zero-points 
by fitting galaxy templates to bright elliptical galaxies in this 
cluster, using BPZ (Berutez 2000) in the Subaru three-band 
photometry combined with HST/ACS images taken in the 

' http://smoka.nao.ac.jp. 
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F435W, F625W, and F850LP bands, for which accurate zero- 
points are well known. Stars were also compared in both Sub- 
aru and the HST/ACS images. Finally, we find consistency of 
our magnitude zero-points to within ±0.05 mag. 

3.2. Weak Lensing Distortion Analysis 

We use the IMC AT package developed by N. Kaiser'' to per- 
form object detection, photometry and shape measurements, 
following the formalism outlined in Kaiser et al. (1995, KSB). 
Our analysis pipeline is implemented based on the procedures 
described in Erben et al. (2001) and on verification tests with 
STEPl and STEP2 data of mock ground-based observations 
(Heymans et al. 2006; Massey et al. 2007). For details of 
our implementation of the KSBh- method, see also Umetsu & 
Broadhurst (2008); Umetsu et al. (2009). 

Object Detection 

Objects are first detected as local peaks in the co-added 
mosaic image in /?c by using the IMCAT hierarchical peak- 
finding algorithm hfindpeaks which for each object yields ob- 
ject parameters such as a peak position, x, an estimate of the 
object size, rg (where rg is the Gaussian scale length of the ob- 
ject as given by hfindpeaks; see KSB for details), the signifi- 
cance of the peak detection, v. The local sky level and its gra- 
dient are measured around each object using the IMCAT get- 
sky routine. Total fluxes and half-light radii, r/,, are then mea- 
sured on sky-subtracted images using the IMCAT apphot rou- 
tine. For details of our IMCAT-based photometry, see §4.2.1 
of Umetsu et al. (2009). We removed from our detection cat- 
alog extremely small objects with rg < \ pixel, objects with 
low detection significance, v < IQ, objects with large raw el- 
lipticities, \e\ > 0.5 (see Clowe et al. 2000; Heymans et al. 
2006; Massey et al. 2007) where |e| = {a^-b^)/{a^ + b^) with 
axis-ratio < fc/a < 1 for images with elliptical isophotes, 
noisy detections with unphysical negative fluxes, and objects 
containing more than 10 bad pixels, nbad > 10. This selec- 
tion procedure yields an object catalog with N = 82431 (93.2 
arcmin"2). 

Weak Lensing Distortion Measurements 

In the KSB algorithm we measure the complex image el- 
lipticity ea = {fill - 222, 612} / (6i 1 + G22) from the weighted 
quadrupole moments of the surface brightness of individual 
galaxies, 

Qaf) = j Sew{e)6o,6pi{e) («,/? = 1,2) (13) 

where I{d) is the surface brightness distribution of an object, 
and W{6) is a Gaussian window function matched to the size 
of the object. In equation (13) the object centroid is chosen as 
the coordinate origin, and the maximum radius of integration 
is chosen to be ^max = 4rg from the centroid. In our practi- 
cal implementation of equation (13), we iteratively refine the 
Gaussian-weighted centroid of each object to accurately mea- 
sure the object shapes. An initial guess for the centroid is pro- 
vided with the IMCAT hfindpeaks routine (see 3.2.0). During 
the process objects with offsets in each iteration larger than 3 
pixels are removed (Oguri et al. 2009). 

The next step is to correct observed image ellipticities for 
the point-spread function (PSF) anisotropy using a sample of 
stars as references: 

^a=^a-Cf«^ (14) 
* http://www.ifa.hawaii.edurkaiser/iincat. 



where Psm is the smear polarizability tensor which is close to 
diagonal, and ^* = iP*m)'^0S* is the stellar anisotropy kernel. 
We select bright (20 ^ ^ 22), unsaturated stellar objects 
identified in a branch of the r/, versus /?c diagram to mea- 
sure ^* . In order to obtain a smooth map of ^* which is 
used in equation (14), we divided the co-added mosaic im- 
age of ~ lOK X 7.7K pixels into 4x3 rectangular blocks, 
where the block length is based on the coherent scale of PSF 
anisotropy patterns (see, e.g., Umetsu & Broadhurst 2008; 
Umetsu et al. 2009). In this way the PSF anisotropy in in- 
dividual blocks can be well described by fairly low-order 
polynomials. We then fitted the q* in each block indepen- 
dently with second-order bi-polynomials, q"(6), in conjunc- 
tion with iterative outlier rejection on each component of the 
residual: Sel = e*^- (P*,,)"'^^* (0). The final stellar sam- 
ple contains 672 stars (i.e., A^* ~ 60 stars per block), or the 
mean surface number density of n* ~ 0.75 arcmin"-^, with 
ly ^ 800. We note that the mean stellar ellipticity before cor- 
rection is e~i* ~ -3.5 X lO"-' and e"2* ~ -5.0 x 10"^ over the 
data field, while the residual e* after correction is reduced to 
g*res^(_0_07±i.i8)x 10-4 and e^''«=' = (-i-2.10±1.58)x 10-4. 
TTie mean offset from the null expectation is reduced down 
to = (2.1 ± 1.6) X 10-4, which is almost two-orders of 

magnitude smaller than the weak lensing signal in cluster out- 
skirts. We show in Figure 2 the quadrupole PSF anisotropy 
field as measured from stellar eUipticities before and after 
the anisotropic PSF correction. Figure 3 shows the distribu- 
tions of stellar ellipticity components before and after the PSF 
anisotropy correction. From the rest of the object catalog, we 
select objects with > + uhairl) as a /?c-selected weak 
lensing galaxy sample, where r,* is the median value of stellar 
half-light radii rjj, corresponding to half the median width of 
the circularized PSF over the data field, cr(rp is the rms dis- 
persion of r^, and is conservatively taken to be 1.5 in the 
present weak lensing analysis. 

We then need to correct image ellipticities for the isotropic 
smearing effect caused by atmospheric seeing and the window 
function used for the shape measurements. The pre-seeing 
reduced shear ga (a = 1,2) can be estimated from 

ga = iP~')afie'^ (15) 

with the pre-seeing shear polarizability tensor defined as 
Hoekstraetal. (1998), 

ng _ nsh [nsmfnsm*\-lnsh*'\ ^ nsh psm "1-^ -I 
nxft-^cx^-y^ ) ^ Ja/3 ~'^«/3~'^«/3tr[psm*] ^^^> 

with P'* being the shear polarizability tensor. In the sec- 
ond equality we have used a trace approximation to the stellar 
shape tensors, P^''* andP™". To apply equation (15) the quan- 
tity tr[P*''*]/tr[P*'"*] must be known for each of the galaxies 
with different size scales. Following Hoekstra et al. (1998), 
we recompute the stellar shapes P'** and P™* in a range of 
filter scales r^ spanning that of the galaxy sizes. At each filter 
scale rg, the median (tr[P'''*]/tr[P""*]) over the stellar sam- 
ple is calculated, and used in equation (16) as an estimate 
of tr[P'**]/tr[P""*]. Further, we adopt a scalar correction 
scheme, namely 

(^g)a/5 = ^tr[Pj],5„^=P^5„^ (17) 

(Erben et al. 2001 ; Hoekstra et al. 1998; Umetsu & Broadhurst 
2008; Umetsu et al. 2009). 
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Following the prescription in Umetsu & Broadhurst (2008) 
and Umetsu et al. (2009), we compute for each object the vari- 
ance of g = gi+ ig2 from neighbors identified in the ob- 
ject size (rg) and magnitude (Re) plane. We take = 30. The 
dispersion Ug is used as an rms error of the shear estimate for 
individual galaxies. With the dispersion Ug we define for each 
object the statistical weight Ug by 



where ai is the softening constant variance (e.g., Hamana 
et al. 2003). We choose ctg = 0.4, which is a typical value 
of the mean rms over the background sample (see, e.g., 
Umetsu & Broadhurst 2008; Umetsu et al. 2009). The case 
with Qf, = corresponds to an inverse-variance weighting. On 
the other hand, the limit Qg ^ CTg yields a uniform weighting. 
Figure 4 shows the mean statistical weight (mj) as a function 
of object size rg (left) and of object magnitude (right) for 
the /?c-selected weak-lensing galaxy sample, where in each 
panel the mean weight is normalized to unity in the first bin. 

Shear Calibration 

In practical observations, the measurement of f| is quite 
noisy for an individual faint galaxy. Further, Pg depends non- 
linearly on galaxy shape moments (see, e.g., Bartelmann & 
Schneider (2001) for an explicit expression for the shape ten- 
sors), so that this nonlinear error propagation may lead to 
a systematic bias in weak lensing distortion measurements 
(T. Hamana, in private communication). Indeed, it has been 
found that smoothing noisy Pg estimates does not necessarily 
improve the shear estimate but can lead to a systematic un- 
derestimate of the distortion amplitude by ~ 10%-15% (e.g., 
Erben et al. 2001; Heymans et al. 2006; Massey et al. 2007). 
In order to improve the precision in shear recovery, we have 
adopted the following shear calibration strategy: First, we se- 
lect as a sample of shear calibrators those galaxies with a de- 
tection significance greater than a certain limit i^^ and with 
a positive raw P| value. Note that the shear calibrator sample 
is a subset of the target galaxy sample (see §3.2.0 and §3.2.0). 
Second, we divide the calibrator size (rg) and magnitude (Rc) 
plane into a grid of 10 x 3 cells^ each containing approxi- 
mately equal numbers of calibrators, and compute a median 
value of Pg at each cell in order to improve the signal-to-noise 
ratio. Finally, each object in the target sample is matched to 
the nearest point on the (rg,Rc) calibration grid to obtain the 
filtered Pg measurement. In this study, we take i^c = 30 as the 

significance threshold for the calibrator sample.^ Finally, we 
use the estimator ga = e'^j (Pg) for the reduced shear. The 
final 7?c-selected galaxy sample contains A' = 38758 objects, 
or the mean surface number density of Ug = 43.7 galaxies 
arcmin"^. The mean variance over the galaxy sample is ob- 
tained as ^ ~ 0.20, or Wg = (^)'/2 ~ 0.46. 

We have tested our weak lensing analysis pipeline includ- 
ing the galaxy selection procedure using simulated Subaru 
Suprime-Cam images of the STEP2 project (Massey et al. 
2007). We find that we can recover the weak lensing signal 
with good precision: typically, m ~ -5% of the shear calibra- 
tion bias, and c ~ 10"^ of the residual shear offset which is 

' We note that Pg is a strong function of the object size rg (see, e.g., Figure 
3 of Umetsu & Broadhurst (2008)), and only weakly depends on the object 
magnitude through 

* The significance threshold for the target sample is 10 (see §3.2.0) 



about one-order of magnitude smaller than the weak-lensing 
signal in cluster outskirts (\g\ ^ 10"^). This level of calibra- 
tion bias is subdominant compared to the statistical uncer- 
tainty (^ 15%) due to the intrinsic scatter in galaxy shapes, 
and the degree of bias depends on the quality of the PSF in 
a complex manner. Therefore, we do not apply shear cal- 
ibration corrections to our data. Rather we emphasize that 
the most critical source of systematic uncertainty in cluster 
weak lensing is dilution of the distortion signal due to the in- 
clusion of unlensed foreground and cluster member galaxies, 
which can lead to an underestimation of the true signal for 
R ^ 400kpc/!"' by a factor of 2-5 (see Figure 1 of Broad- 
hurst et al. 2005b). 

3.3. Galaxy Sample Selection from the Color-Color Diagram 

It is crucial in the weak lensing analysis to make a secure se- 
lection of background galaxies in order to minimize contami- 
nation by unlensed cluster and foreground galaxies; otherwise 
dilution of the distortion signal arises from the inclusion of 
unlensed galaxies, particularly at small radius where the clus- 
ter is relatively dense (Broadhurst et al. 2005b; Medezinski 
et al. 2007; Umetsu & Broadhurst 2008). This dilution effect 
is simply to reduce the strength of the lensing signal when 
averaged over a local ensemble of galaxies, in proportion to 
the fraction of unlensed galaxies whose orientations are ran- 
domly distributed, thus diluting the lensing signal relative to 
the true background level, derived from the uncontaminated 
background level we will establish below. 

To separate unlensed galaxies from the background and 
hence minimize the weak lensing dilution, we follow 
the background selection method recently developed by 
Medezinski et al. (2009), which relies on empirical correla- 
tions in color-color (CC) space derived from the deep Subaru 
photometry, by reference to the deep photometric redshift sur- 
vey in the COSMOS field (see §3.4). For C10024H-1654, we 
have a wide wavelength coverage (BR^z') of Subaru Suprime- 
Cam. When defining color samples, we require that objects 
are detected in all three Subaru bands. Further, we hmit the 
data to z' ^ 25.5 AB mag in the reddest band available for 
the cluster. Beyond this Umit incompleteness creeps into the 
bluer bands, complicating color measurements, in particular 
of red galaxies. Our CC-selection criteria yielded a total of 
A' = 8676, 1655, and 5004 galaxies for the red, green, and 
blue samples, respectively, usable for our weak lensing distor- 
tion analysis; these correspond to mean surface number densi- 
ties of fig = 10.9,2.1, and 6.3 galaxies arcmin"^, respectively. 
In Table 3, we list the magnitude limits, the total number of 
galaxies (A'), the mean surface number density (fig), and the 
mean rms error for the galaxy shear estimate (ffg), for our 
color samples. The resulting color boundaries for respective 
galaxy samples are shown in Figure 5. We emphasize that this 
model-independent empirical method allows us to clearly dis- 
tinguish these distinct blue and red populations (see Figures 5 
and 6), by reference to the well-calibrated COSMOS photom- 
etry, as well as to separate unlensed foreground/cluster galax- 
ies from the background (see also §3.4 and Figure 7). 

Cluster Galaxies 

Following the prescription by Medezinski et al. (2009) we 
construct a CC diagram and first identify in this space where 
the cluster lies by virtue of the concentration of cluster mem- 
bers. Figure 6 (left panel) shows in CC space the distribu- 
tion of mean projected distances (9) from the cluster center 
for all galaxies in the C10024-H1654 field. As demonstrated 
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by Medezinski et al. (2009), this diagram clearly reveals the 
presence of strong clustering of galaxies with lower mean ra- 
dius, confined in a distinct and relatively well-defined region 
of CC space. This small region corresponds to an overden- 
sity of galaxies in CC space comprising the red sequence of 
the cluster and a blue trail of later type cluster members, as 
clearly demonstrated in the right panel of Figure 6. In Figures 
5 (green points) and 6 (left panel) we mark this overdense 
region in CC space to safely encompass the region signifi- 
cantly dominated by the cluster. We term the above sample, 
which embraces all cluster member galaxies, the green sam- 
ple, as distinct from well separated redder and bluer galaxies 
identified in this CC space (see the right panel of Figure 6). 
Naturally, a certain fraction of background galaxies must be 
also expected in this region of CC space, where the proportion 
of these galaxies can be estimated by comparing the strength 
of their weak lensing signal with that of the reference back- 
ground samples (Medezinski et al. 2007, 2009). Figure 10 
(crosses) demonstrates that the level of the tangential distor- 
tion for the green sample is consistent with zero to the out- 
skirts of the cluster, indicating that the proportion of back- 
ground galaxies in this sample is small compared with cluster 
members. 

Red Background Galaxies 

To define the foreground and background populations, we 
utilize the combination of the strength of the weak lensing 
signal and the number density distribution of galaxies in CC 
space. With the BR^z' photometry we can improve upon the 
simple color-magnitude selection previously performed in our 
weak lensing analyses of z ~ 0.2 clusters (Medezinski et al. 
2007; Umetsu & Broadhurst 2008; Umetsu et al. 2009), where 
we had defined a "red" population, which comprises mainly 
objects lying redward of the E/SO color-magnitude sequence. 
This can be done by properly identifying and selecting in CC 
space the reddest population dominated by an obvious over- 
density and a red trail (see the right panel of Figure 6). For 
this red sample we define a conservative diagonal boundary 
relative to the green sample (see Figure 5, red points), to 
safely avoid contamination by cluster members (§3.3.0) and 
also foreground galaxies (see § 3.3.0). To do this, we measure 
the average distortion strength as a function of the distance 
from the cluster sequence in CC space, and take the limit to 
where no evidence of dilution of the weak lensing signal is 
visible. We also define a color limit at R^-zf — 0.35 that sep- 
arates what appears to be a distinct density maximum of very 
blue objects. The boundaries of the red sample as defined 
above are marked on Figure 5 (red points), and can be seen 
to lie well away from the green cluster sample. For this red 
sample we show below a clearly rising weak-lensing signal 
all the way to the smallest radius accessible (Figure 10, trian- 
gles), with no sign of a central turnover which would indicate 
the presence of unlensed cluster members. 

Blue Background and Foreground Galaxies 

Special care must be taken in defining blue background 
galaxies, as objects lying bluer than the E/SO sequence of a 
cluster can comprise blue cluster members, foreground ob- 
jects, and background blue galaxies. This is of particular con- 
cern where only one color (i.e., 2 bands) is available. We have 
shown in our earlier work (Broadhurst et al. 2005b; Medezin- 
ski et al. 2007) that this can lead to a dilution of the weak 
lensing signal relative to the red background galaxies due to 



unlensed foreground and cluster galaxies, the relative propor- 
tion of which will depend on the cluster redshift. Encourag- 
ingly, it has been demonstrated by Medezinski et al. (2009) 
that the foreground unlensed population is well defined in CC 
space as a clear overdensity (Figure 5, magenta points; Figure 
6, right panel) and we can therefore simply exclude these ob- 
jects in this region of CC space from our analysis by setting 
the appropriate boundaries relative to this overdensity, found 
to be where the weak lensing signal starts showing dilution 
by these foreground galaxies. The bluer galaxy overdensity 
in CC space, seen in Figures 5 and 6, is also unclustered (Fig- 
ure 5, blue points) and mainly concentrated in one obvious 
cloud in CC space. This blue cloud has a continuously rising 
weak lensing signal (Figure 10, circles), towards the center of 
the cluster, with an amplitude which is consistent with the red 
background population defined and hence we can safely con- 
clude that these objects lie in the background with negligible 
cluster or foreground contamination, which would otherwise 
drag down the central weak lensing signal. The boundaries 
of this blue background sample are plotted in Figure 5 (blue 
points) which we extend to include object lying outside the 
main blue cloud but well away from the foreground and clus- 
ter populations defined above. 

3.4. Depth Estimation 

An estimate of the background depth is required when con- 
verting the observed lensing signal into physical mass units, 
because the lensing signal depends on the source redshifts 
through the distance ratio /3(zs) = Djs/Ds. 

Since we cannot derive complete samples of reliable pho- 
tometric redshifts from our limited three-band (BR^z') images 
of C10024H-1654, we instead make use of deep field photome- 
try covering a wider range of passbands, sufficient for photo- 
metric redshift estimation of faint field redshift distributions, 
appropriate for samples with the same color-color/magnitude 
limits as our red and blue populations. The 30-band COS- 
MOS photometry (Ilbert et al. 2009) is very suited for our 
purposes, consisting of deep optical and near-infrared pho- 
tometry over a wide field, producing reliable photometric red- 
shifts for the majority of field galaxies to faint Umiting magni- 
tudes: m < 25 AB mag in the Subaru i' band. The public 30- 
band COSMOS photometric catalog contains about 380000 
objects over 2deg^ covering Subaru BVg'r'i'z' photometry. 
The photometric zero-point offsets given in Table 13 of Ca- 
pak et al. (2007) were applied to the COSMOS catalog. From 
this we select ~ 3 x 10^ galaxies with reliable photometric 
redshifts as a COSMOS galaxy sample. Since the COSMOS 
photometry does not cover the Subaru R^ band, we need to 
estimate Rc-bmd magnitudes for this COSMOS galaxy sam- 
ple. We use a new version of the HyperZ template fitting 
code (New-HyperZ ver.ll, Roser Pello, private communica- 
tion; Bolzonella et al. 2000) to obtain for each galaxy the best- 
fitting spectral template, from which the R^ magnitude is de- 
rived with the transmission curve of the Subaru /?c-band filter. 
Note, since Subaru BVg'r'i'z' magnitudes are available for all 
these galaxies, this Rc estimation can be regarded as an in- 
terpolation. Therefore, the /?c magnitudes obtained with this 
method will be sufficiently accurate for our purpose, even if 
photometric redshifts derived by HyperZ (which will not be 
used for our analysis) suffer from catastrophic errors. 

To assess the effective redshift depth for our blue and red 
background populations, we apply for each sample our color- 
color/magnitude selection to the COSMOS multiband pho- 
tometry, and obtain the redshift distribution Niz) of the back- 
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ground population with the same color and magnitude cuts. 
The resulting photometric -redshift distributions N{z) of the 
CC/magnitude-selected red, green, and blue samples in the 
COSMOS field are displayed in Figure 7 (see also Medezin- 
ski et al. 2009). We then calculate moments of the redshift 
distribution of the distance ratio /3(zj) = D^s/Ds for each back- 
ground population as 



(/?"> = 



/ dzN(z),l3"(z) 
J dzN(z) 



(19) 



The first moment (/?) represents the mean lensing depth in the 
weak lensing limit (k, \'y\ ^ 1, g ~ 7), where the relationship 
between the surface mass density and the lensing observables 
is linear (see §2). In this case, one can safely assume that all 
background galaxies lie in a single source plane at redshift 
Zs,^ corresponding to the mean depth defined as {(3) = 
Pizs,0)- Table 3 lists for respective color samples the mean 
source redshift Zs, the effective single-plane redshift Zs,^, and 
the first and second moments of the distance ratio. From this, 
we find the blue sample to be deeper than the red by a factor 
of (/3(Blue)) /(/3(Red)) = 1.21 ±0.08, which is consistent with 
the corresponding ratio, (g+(Blue))/(g+(Red)) = 1.17 ±0.20 
(see Figure 10), of the mean tangential distortion averaged 
over the full radial extent of the cluster. 

In general, a wide spread of the redshift distribution of 
background galaxies, in conjunction with the single plane ap- 
proximation, may lead to an overestimate of the gravitational 
shear in the nonlinear regime (Hoekstra et al. 2000). To the 
first order of k, this bias in the observed reduced shear is writ- 
ten as (Hoekstra et al. 2000; Seitz & Schneider 1997) 
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(20) 



where g is the reduced shear from the single source plane 
assumption, namely, g = gds.fs) = "/(Zs.i3)/[l - nds.js)], and 
K = K,izs,i3)- With the 30-band COSMOS photometry, the level 
of bias is estimated as A.g/g w 0.020k, 0.016k, and 0.021k 
for the red, blue, and blue+red background samples, respec- 
tively. 

By virtue of the great depth of the Subaru imaging and of 
the moderately low redshift of the cluster, for C10024+I654, 
this effect turns out to be quite negligible at all radii in the sub- 
critical regime (k < 1). Finally, taking into account the photo- 
metric zero-point errors in our Subaru photometry (§3.1) and 
the photometric redshift errors of individual COSMOS galax- 
ies, we estimate the uncertainty in the mean depth (/3) to be 
~ 4% for the red galaxies, and ~ 6% for the blue galaxies; for 
the composite blue-nred background sample with a blue-to-red 
ratio of 0.58 (see Table 3), it is about 5%. 

4. CLUSTER WEAK LENSING ANALYSIS 
4.1. Two-Dimensional Mass Map 

Weak lensing measurements of the gravitational shear field 
can be used to reconstruct the underlying projected mass den- 
sity field, S„(0). In the present study, we use the dilution- 
free, B7?cz'-selected blue+red background sample (§3.3) both 
for the 2D mass reconstruction and the lens profile measure- 
ments. We follow the prescription described in Umetsu et al. 
(2009, see §4.4) to derive the projected mass distribution of 
the cluster from Subaru distortion data. 

Figure 8 displays the spin-2 reduced-shear field, g{0), ob- 
tained from the blue+red background sample, where for vi- 
sualization purposes g{d) is smoothed with a Gaussian with 



1.41' FWHM. In the left panel of Figure 9 we show the 
reconstructed two-dimensional map of lensing convergence 
K{d) = E„(0)/Scrit in the central 26' x 22' region. A promi- 
nent mass peak is visible in the cluster center, around which 
the lensing distortion pattern is clearly tangential (Figure 8). 
This first maximum in the k map is detected at a significance 
level of 16(7, and coincides well with the optical cluster cen- 
ter within the statistical uncertainty: AR.A. = 4.5" ± 4.6", 
ADecl. = 0.0" ±5.6", where AR.A. and ADecl. are right- 
ascension and decUnation offsets, respectively, from the cen- 
ter of the central bright elliptical galaxy, or the galaxy 374 in 
the spectroscopic catalog of Czoske et al. (2002). Our mass 
peak is also in spatial agreement with the X-ray emission peak 
revealed by Chandra ACIS-S observations of Ota et al. (2004, 
Gl peak), R.A. ~ 00:26:36.0, Deck ~ +17:09:45.9 (J2000.0), 
which is offset to the northeast by ~ 5" from the optical cen- 
ter, and is close to the galaxy 380 in Czoske et al. (2002). 

Also compared in Figure 9 are member galaxy distributions 
in C10024+1654, Gaussian smoothed to the same resolution 
of FWHM = 1.41'. The middle and right panels display the 
/?c— band number and luminosity density maps, respectively, 
of green cluster galaxies (see Table 3). Overall, mass and light 
are similarly distributed in the cluster: The cluster is fairly 
centrally concentrated in projection, and contains an extended 
substructure located about 3' (~ 670 kpc h~^) northwest of the 
cluster center, as previously found by CFHTAVHT spectro- 
scopic observations of Czoske et al. (2002) and HST/WFPC2 
and STIS observations of Kneib et al. (2003). This NW 
clump of galaxies is associated with the primary cluster in 
redshift space (peak A at z = 0.395 in Czoske et al. (2002), 
and the density peak of the NW galaxy clump is located at 
AR.A. ~ -2.29', ADecl. ~ 2.48' in our galaxy number den- 
sity map. 

4.2. Cluster Center Position 

In order to obtain meaningful radial profiles one must 
carefully define the center of the cluster For this purpose 
we rely on the improved strong-lens model of Zitrin et al. 
(2009b), constructed using deep HST/ACS and NIC3 im- 
ages. Their mass model accurately reproduces the well- 
known, spectroscopically-confirmed 5-image system of a 
source galaxy at z = 1 .675 (Broadhurst et al. 2000), and con- 
firms the tentative 2-image system identified by Broadhurst 
et al. (2000), finding an additional third image associated 
with this source. In addition, 9 new multiple-image sys- 
tems were identified by their improved mass model, bringing 
the total known for this cluster to 33 multiply-lensed images 
spread fairly evenly over the central region, 6 ^AW . The 
mass model of Zitrin et al. (2009b) reveals a fairly round- 
shaped radial critical curve with radius ~ 10" (at z = 1.675), 
providing a reasonably well defined center (AR.A. ~ -1.5", 
ADecl. ^ -2.5"), which is slightly offset, but located fairly 
close to, the optical center, and was used as the center of mass 
in the radial mass-profile analysis of Zitrin et al. (2009b). In 
this study, we will define the center of mass in a more quanti- 
tative manner as the peak position of the smooth dark-matter 
component of the Zitrin et al. mass model. The resulting cen- 
ter of mass (dark-matter center, hereafter) is at offset position 
AR.A. ~ -2.32", ADecl. ~ -1.44", consistent with the geo- 
metric center of the inner critical curve within 1.3". We note 
that the mass peak in our Subaru k map is in spatial agreement 
with this dark-matter center. 

The central mass distribution of C10024+1654 has been ex- 
amined by other authors using high-resolution HST observa- 
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tions. Kneib et al. (2003) determined the center of mass to be 
at AR.A. ~ -2.3", ADecl. ~ -4.2", from their WFPC2/STIS 
weak lensing measurements and strong lensing constraints 
by the 5-image system (z = 1 .675), and their center of mass 
is close to our dark-matter center ~ 2.7"). Jee et al. 
(2007) derived a high-resolution mass map from their lens- 
ing analysis of deep 6-band ACS images, incorporating as 
strong-lensing constraints the 5-image system and two addi- 
tional multiple-system candidates (Objects B1-B2 and C1-C2, 
in their notation).'' Their mass peak is in good agreement with 
central bright elliptical galaxies, and close to the galaxy 374 
in Czoske et al. (2002). Jee et al. (2007) chose the geomet- 
ric center of the ringlike dark-matter structure (at 6 ^ 75"), 
revealed by their non-parametric mass reconstruction, as the 
cluster center for their radial mass-profile analysis. Their ring 
center is located at AR.A. ~ 3.3", ADecl. ~ -7.6", and is 
about 6.5" offset from our dark-matter center. This offset is 
about 9% of the dark-matter ring radius (9 ^ 75"); this level 
of discrepancy can be reconciled by noting that the ringlike 
structure revealed by Jee et al. (2007) is diffuse and not per- 
fectly round in shape. 

4.3. Lens Distortion Profile 

The spin-2 shape distortion of an object due to gravita- 
tional lensing is described by the complex reduced shear, 
g = gi +ig2 (see equation [10]), which is coordinate depen- 
dent. For a given reference point on the sky, one can instead 
form coordinate-independent quantities, the tangential distor- 
tion g+ and the 45° rotated component, from linear combina- 
tions of the distortion coefficients and g2 as 

g+ = -(^iCOs2(/)+^2sin2(/)), = -(^2cos20-gi sin2(/)), 

(21) 

where is the position angle of an object with respect to 
the reference position, and the uncertainty in the g+ and gx 
measurement is a+ = ax = Cg/v^ = cr in terms of the rms 
error ag for the complex shear measurement. Following the 
ACS strong lensing analysis of Zitrin et al. (2009b), we take 
their dark-matter center as the cluster center of mass for our 
radial profile analysis (see §4.2). To improve the statistical 
significance of the distortion measurement, we calculate the 
weighted average of g+ and g 

.J2iUg,i8+,i 



{gAOm))-- 



(22) 



(23) 



where the index i runs over all of the objects located within 
the mth annulus, and Kg,, is the statistical weight (see equation 
[18]) for the ith object, and 6m is the weighted center of the 
mth radial bin. 
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(24) 



with 9i being the offset vector of the ;th galaxy position from 
the cluster center (§4.2). In our practical analysis, we use the 
continuous limit of equation (24): See Appendix A of Umetsu 
& Broadhurst (2008). 

In Figure 10 we compare azimuthally-averaged radial pro- 
files of the tangential distortion g+ (E mode) and the 45° ro- 
tated (x) component gx (B mode) as measured for our red, 

' Regarding the validity of the lensing hypothesis of the C1-C2 system, 
see discussions in Zitrin et al. (2009b). 



blue, and green galaxy samples (Table 3). The error bars rep- 
resent 68.3% confidence limits, (t+, estimated by bootstrap 
resamphng of the original data set. The red and blue popula- 
tions show a very similar form of the radial distortion profile 
which declines smoothly from the cluster center, remaining 
positive to the limit of our data, 6'max =16'. The mean distor- 
tion amplitude of the blue population is consistent with, but 
slightly higher than, that of the red population, which is re- 
lated to the greater depth of the blue population relative to the 
red (§3.3 and 3.4). This smooth overall trend suggests that 
the NW substructure identified in Figure 9 has only a minor 
effect on the overall profile, as found in the WFPC2 weak- 
lensing analysis of Kneib et al. (2003). On the other hand, 
the tangential distortion profile for the green galaxies is con- 
sistent with a nuU signal at aU radii, while this population is 
strongly clustered at small radius (Figures 6 and 9), indicat- 
ing that the green galaxies mostly consists of cluster member 
galaxies. This convincingly demonstrates the power of our 
color selection method. 

Now we assess the tangential distortion profile from the 
blue-Hred background sample (§3.3) in order to examine the 
form of the underlying cluster mass profile and to characterize 
cluster mass properties. In the weak lensing limit (k, |7| ^ 1), 
the azimuthally averaged tangential distortion profile {g+iO)) 
(equation [22]) is related to the projected mass density profile 
(e.g., Bartelmann & Schneider 2001) as 



{g^ie))^{^Ae)) = Ri<e)-{m), 



(25) 



where (• • •) denotes the azimuthal average, and k(< 9) is the 
mean convergence within a circular aperture of radius 9 de- 
fined as K« 9) = (7r6i2)-i J^Q'^^g(fe' k(0'). Note that the sec- 
ond equality in equation (25) holds for an arbitrary mass dis- 
tribution. With the assumption of quasi-circular symmetry in 
the projected mass distribution, one can express the tangential 
distortion as 

in the non-linear but sub-critical {AetA{d) > 0) regime. Fig- 
ure 11 shows the tangential and 45°-rotated distortion pro- 
files for the blue-nred background sample. Here the pres- 
ence of B modes can be used to check for systematic errors. 
The observed ii-mode signal is significant with a total detec- 
tion signal-to-noise ratio (S/N) of ~ 14, remaining positive 
to the limit of our data, 9^^^ = 16', or a projected distance of 
3.6Mpc/z~'. The x -component is consistent with a null sig- 
nal at all radii, indicating the reliability of our weak lensing 
analysis. The significance level of the B-mode detection is 
about 2.8tT, which is about a factor of 5 smaller than £-mode. 

To quantify and characterize the mass distribution of 
C10024H-1654, we compare the measured g+ profile with the 
physically and observationally motivated NFW model (e.g., 
Broadhurst et al. 2005b; Umetsu &. Broadhurst 2008; Broad- 
hurst et al. 2008; Umetsu et al. 2009; Okabe et al. 2009). The 
NFW universal density profile has a two-parameter functional 
form as 



PNFW('") ■■ 



{r/r,){\ + r/rsf' 



(27) 



where ps is a characteristic inner density, and is a character- 
istic inner radius. The logarithmic gradient n = d\n p(r) /d]nr 
of the NFW density profile flattens continuously towards the 
center of mass, with a flatter central slope « = -1 and a steeper 
outer slope (n — )• -3 when r — > oo) than a purely isother- 
mal structure (n = -2). A useful index, the concentration, 
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compares the virial radius, rvir, to of the NFW profile, 
Cvir = ^vir/rv- We Specify the NFW model with the halo virial 
mass Myir and the concentration Cm instead of ps and r^. 
We employ the radial dependence of the NFW lensing pro- 
files, knfw(^) and 7+_nfw(^^), given by Bartelmann (1996) and 
Wright & Brainerd (2000).'" The NFW density profile can 
be further generalized to describe a profile with an arbitrary 
power-law shaped central cusp, p(r) oc r~", and an asymptotic 
outer slope of n = -3 (Zhao 1996; Jing & Suto 2000), 

PgNFw(r) = I , , (28) 

which reduces to the NFW model for a= I. We refer to 
the profile given by equation (28) as the generalized NFW 
(gNFW) profile. It is useful to introduce the radius r_2 at 
which the logarithmic slope n of the density is isothermal, 
i.e., n = -2. For the gNFW profile, r_2 = (2 - a)rs, and 
thus the corresponding concentration parameter reduces to 
c_2 = r^ir/r-i = Cvir/(2-a). We specify the gNFW model with 
the central cusp slope, a, the halo virial mass, Mvir, and the 
concentration, c-2 = Cnr/Q. — a). 

Table 4 summarizes the results of fitting with the NFW 
(gNFW) model, listing the lower and upper radial limits of 
the data used for fitting, the best-fit parameters (Myir, Cyir) and 
their errors (68.3% CL), the minimized value (x^) with 
respect to the degrees of freedom (dof), and the predicted Ein- 
stein radius for a background source at Zs = 1.675, cor- 
responding to the 5-image system. The quoted errors in Ta- 
ble 4 include the uncertainty in the mean depth (/?) of the 
background galaxies (see Table 3). The best-fit NFW model 
is given by My„ = 1.14!° x lO'^M© r' and Cvir = 10.6!|g, 
with x^n/dof = 2.2/8. This model yields an Einstein radius 
of = 'iT'^ii at Zs = 1.675, consistent within the errors with 
the observed location of the Einstein radius ^e = 30" ± 3" 
(z, = 1.675, Zitrin et al. 2009b). 

Assuming a singular isothermal sphere, this Einstein ra- 
dius constraint is readily translated into the equivalent one- 
dimensional velocity dispersion of cTv = 1245 ± 62km s"'. 
From a fit of the truncated, nonsingular isothermal sphere pro- 
file to the strong-lensing mass model of Tyson et al. (1998), 
Shapiro & Iliev (2000) obtained an average velocity disper- 
sion of (jy ~ 1200 km s~' within a sphere of r = 600kpc/!~', 
in close agreement with the value of ctv = 1 150 ± 100 km s"' 
measured by Dressier et al. (1999) from 107 galaxy red- 
shifts. A similar value of cr,, = 1050 ±75 km s"' was de- 
rived by Czoske et al. (2001, 2002) from 193 galaxy redshifts 
within 5' from the cluster center. Diaferio et al. (2005) used 
their caustic method to study the internal velocity structure of 
C10024+1654. This method allows one to identify in a radius- 
redshift phase space diagram the cluster boundaries, which 
serve as a measure of the local escape velocity (see Diafe- 
rio 1999; Lemze et al. 2008). With the spectroscopic galaxy 
catalog of Czoske et al. (2001, 2002) they found an average 
one-dimensional velocity dispersion of ~ 937 km s~' for 

Note that the Bartelmann's formula for the NFW lensing profiles are 
obtained assuming that the projection integral to extend to infinity. Alter- 
natively, a truncated NFW profile can be used to model the lensing pro- 
filesTakada & Jain (2003). We have confirmed that the best-fit NFW param- 
eters obtained using the above two different formulae agree to within 0.5% 
for the case of C10024-H1654 lensing; for detailed discussions, see Baltz et al. 
(2007); Hennawi et al. (2007). 

" For a given source redshift, the Einstein radius (?h is defined as 1 = k(< 
Oa-Zs). For an NFW model, this equation for 9^ can be solved numerically, 
for example, by the Newton-Raphson method. See Appendix B of Umetsu & 
Broadhurst (2008). 



this cluster, which is lower but consistent with the previous 
results. A comparison between their caustic and our lensing 
mass profiles will be given in §6.2. 

4.4. Lensing Convergence Profile 

Here we examine the lensing convergence (k) profile using 
the one-dimensional, non-parametric reconstruction method 
developed by Umetsu & Broadhurst (2008, Appendix A) 
based on the nonlinear extension of aperture mass densito- 
metry which measures the projected mass interior to a given 
radius from distortion data. See also Appendix of Umetsu 
et al. (2009) for details of this reconstruction method. 

We use a variant of the aperture mass densitometry, or the 
so-called C-statistic (Fahlman et al. 1994; Clowe et al. 2000), 
of the form: 



i-iebi/Obif 

= k(< 0)-Kb, 



Obi 



d\n6' {^^{6')) 



(29) 



where k(< 9) is the average convergence interior to radius 
9, and {9bi , 9b2) are the inner and outer radii, respectively, of 
the annular background region in which the mean background 
contribution, ki, = R{9h\ <'d < 9h2), is defined. The substruc- 
ture contribution to k is local, whereas the inversion from the 
observable distortion to k involves a non-local process (§2). 
Consequently the one-dimensional inversion method requires 
a boundary condition specified in terms of the mean back- 
ground convergence k/,. The inner and outer radii of the an- 
nular background region are set to 6h\ = 14' (3.1 Mpc/z"') and 
9h2 = 16' (3.6 Mpc /;"'), respectively, sufficiently larger than 
the cluster virial radius of massive clusters (^2Mpc/i"'), 
so that the weak-lensing approximation (7+(^)) « {g+i9)) is 
valid in the background region. In the nonlinear but subcrit- 
ical regime (i.e., outside the critical curves), {'y+{9)) can be 
expressed in terms of the averaged tangential reduced shear as 
(7+(^)) PS {g+i9))[l - (k(^))] assuming a quasi-circular sym- 
metry in the projected mass distribution. Hence, for a given 
boundary condition Rh, the non-linear equation (29) for Cc(^) 
can be solved by an iterative procedure: 

Cr"(^)«2 / d\n9'{gm)[l-{K^'\0))] 
Je 

2 



\-{9bi/9b2f Je,, 
{K^'\0)) = CeC^'\0) + -nb 



/ d\ne'{g^{e')), (30) 

Jen 

(31) 



where we have introduced a differential operator defined as 
Co = 252 J^^^ '^'^at satisfies Cel = 1, and C^*' represents 
the aperture densitometry in the kth step of the iteration 
{k = 0, 1,2, ...). This iteration is preformed by starting with 
{Kf-^^{9)) = for all radial bins, and repeated until convergence 
is reached at all radial bins. We compute the bin-to-bin er- 
ror covariance matrix C,y = {5n{9i)5K{9j)) for (k,) = {n{9i)) 
with the linear approximation by propagating the rms error 
cr+{0) for the averaged tangential shear measurement {g+(9)). 
Finally, we determine the best value of Rb iteratively in the 
following way: The iterations start with Rh = 0. At each itera- 
tion, we update the value of Rb using the best-fit NFW model 
(i.e., assuming p{r) oc at large radii) for the reconstructed 



10 



Mass Structure of C10024+1654 from a Full Lensing Analysis 



K profile. This iteration is repeated until convergence is ob- 
tained. 

In Figure 12 we show the resulting surface-mass density 
profile reconstructed from the lens-distortion measurements 
of our blue-nred background galaxies registered in deep Sub- 
aru BR^z' images. The error bars are correlated, with neigh- 
boring bins having ^ 10% cross-correlation coefficients at 
inner radii, and a ~ 40% cross-correlation coefficient at the 
outermost radii. The best-fit NFW model to the k profile is 
obtained as Mvk- = 1.08!|];?^ x lO'^Mo/;-' and = 9.1^^^ 
(X^i„/dof =3.6/8, see Table 4), with the resulting Rb value 
of 2.74 X 10"^, fully consistent with the results from the g+ 
profile (§4.3), ensuring the validity of the boundary condition 
for a shear-based mass reconstruction (Umetsu & Broadhurst 

2008) . We note that simply assuming = yields very sim- 
ilar results, Mvir = 1.06!^;^^ x lO'^Mgr' and Cvir = 10.9!^^3° 
with /dof = 3 .4/8, in agreement well within the 1 a uncer- 
tainty. Our reconstructed k, profile is fairly smooth and well 
approximated by a single NFW profile, but exhibits a sUght 
excess ?ii 6 ^ 3', although with large scatters, with respect to 
the best-fit NFW profile. This excess feature roughly coin- 
cides with the projected distance of the NW clump identified 
in Figure 9. We will come back to this in §5.2. 

From the NFW fit to the k profile the statistical uncer- 
tainty on Mvii- is about 20%, comparable to that from the g+ 
profile, while the constraint on Cvii is rather weak, allowing 
a wide range of the concentration parameter: 5 ^ Cyk ^ 25 
(68.3% CL). This reflects the fact that the reconstruction er- 
ror tends to increase towards the central region as a result of 
inward error propagation. Consequently, the Einstein radius 
is poorly constrained from the one-dimensional mass recon- 
struction: = 33"!!2o Zs = 1.675. Nevertheless, this one- 
dimensional, non-parametric inversion method allows us to 
derive a n profile with a full covariance matrix from weak- 
lensing distortion data alone, which can be readily compared 
and combined with inner strong lensing data to provide a full 
mass profile for the entire cluster (see §5). Furthermore, such 
non-parametric mass profiles are useful when comparing the 
total matter distribution with cluster properties obtained from 
other wavelengths (Lemze et al. 2008, 2009; Lapi & Cavaliere 

2009) . 

4.5. Lensing Depletion Profile 

Lensing magnification, iJL{d), influences the observed sur- 
face density of background sources, expanding the area of 
sky, and enhancing the observed flux of background sources 
(Broadhurst et al. 1995; Umetsu & Broadhurst 2008). The for- 
mer effect reduces the effective observing area in the source 
plane, decreasing the number of background sources per solid 
angle; on the other hand, the latter effect amplifies the flux of 
background sources, increasing the number of sources above 
the limiting flux. 

For the number counts to measure magnification, we use 
our red background sample based on the SExtractor photom- 
etry (§3.1). For these the intrinsic count slope s = d\ogNi< 
m)/dm at faint magnitudes is relatively flat, ^ 0.1, so that a 
net count depletion results (Broadhurst et al. 2005b; Umetsu 
& Broadhurst 2008; Broadhurst et al. 2008). For depletion 
analysis, we have a total of 18561 red galaxies down to a hm- 
iting magnitude of = 25.5 AB mag (see §3.3), where the 
sample size is about twice as large as that for distortion anal- 
ysis; the smaller sample for the distortion analysis is due to 
the fact that it requires the galaxies used are weU resolved to 



make reliable shape measurements (Broadhurst et al. 2005b; 
Umetsu & Broadhurst 2008). 

The number counts for a given magnitude cutoff rucau 
approximated locally as a power-law cut with slope, s = 
t/logjoA^(< m)/dm, are modified in the presence of lensing as 
A^(< nicut) ~ A^o(< ?Mcut)M^'^'"' (Broadhurst et al. 1995), where 
A'o(< nicut) is the unlensed counts. Thanks to the large field 
of view of Subaru/Suprime-Cam, the normalization and slope 
of the unlensed counts for our red galaxy sample are reliably 
estimated as «o = 20.2 ± 0.4arcmin"^ and s = 0.12 ± 0.06, re- 
spectively, from the outer region of the cluster, 11' ^ 9 ^ 16'. 
The slope is less than the lensing invariant slope, s = 0.4, so a 
net deficit of background galaxies is expected. 

The count-in-cell statistic is measured from the flux-Umited 
red background sample on a regular grid of 70 x 52 equal- 
area cells covering a field of 35' x 26' (Umetsu & Broadhurst 
2008). We then calculate the radial profile of the red galaxy 
counts by azimuthally averaging the count-in-cell distribu- 
tion, where each cell is weighted by the fraction of its area 
lying within the respective annular bins (Umetsu & Broad- 
hurst 2008) and a tail of > 3a cells is excluded in each an- 
nulus to remove inherent small-scale clustering of the back- 
ground (Broadhurst et al. 2008). Here the masking effect due 
to bright cluster galaxies and bright foreground objects has 
been properly taken into account and corrected for, following 
the prescription of Umetsu & Broadhurst (2008, see §4.2). 
We conservatively account for the masking of observed sky 
by excluding a generous area nab around each masking ob- 
ject, where a and b are defined as fmask = 3 times the ma- 
jor (A_IMAGE) and minor axes (B_IMAGE) computed from 
SExtractor, corresponding roughly to the isophotal detection 
limit (see Umetsu & Broadhurst 2008). We calculate the cor- 
rection factor for this masking effect as a function of radius 
from the cluster center, and renormalize the number density 
of each radius accordingly. The masking area is estimated as 
about 3% at large radii, and found to increase up to ~ 25% of 
the sky close to the center, R ^ 170kpc/!~'. Note that if we 
use the masking factor i^mask of 2 or 4, instead of 3, the re- 
sults shown below remain almost unchanged (see for details, 
Umetsu & Broadhurst 2008, §5.5.3). 

Figure 13 shows the resulting count-depletion profile de- 
rived from the red background sample based on the SExtractor 
photometry. The error bars include not only the Poisson con- 
tribution but also the variance due to variations of the counts 
along the azimuthal direction, i.e., contributions from the in- 
trinsic clustering of red galaxies and departure from circu- 
lar symmetry (similar to the second term of equation (42) of 
Umetsu & Broadhurst 2008). A strong depletion of the red 
galaxy counts is shown in the central, high-density region of 
the cluster, and clearly detected out to a few arcminutes from 
the cluster center. The statistical significance of the detec- 
tion of the depletion signal is about 8(t. The detection sig- 
nificance of the distortion derived from the blue-nred back- 
ground sample (14(t, see §4.3) is better than the counts, so 
that here we use our depletion measurements only as a consis- 
tency check. The magnification measurements with (squares) 
and without (triangles) the masking correction are roughly 
consistent with each other. To test the consistency between 
our distortion and depletion measurements, we calculate the 
depletion of the counts, n(9) = noM'^'^*~\ expected for the best- 
fitting NFW profile derived from the distortion measurements 
(Figure 11), normalized to the observed density no- A slight 
dip at = 2' - 3' in the depletion profile corresponds to the 
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contribution of the NW clump, which is also seen in the Sub- 
aru distortion data. This comparison shows clear consistency 
between two independent lensing observables with different 
systematics, which strongly supports the reliability and accu- 
racy of our weak lensing analysis. The count depletion of red 
galaxies is seen in all our earlier work (A1689, A1703, A370, 
RXJ1347-11), and the clear result (Scr detection) found here 
strengthens the use of this information when testing for con- 
sistency with weak distortions (e.g., Broadhurst et al. 2005b; 
Umetsu & Broadhurst 2008; Broadhurst et al. 2008). 

5. COMBINING STRONG AND WEAK LENSING DATA 

The Subaru data allow the weak lensing profiles to be ac- 
curately measured in several independent radial bins in the 
subcritical regime (0 > 0^). Here we examine the form of 
the projected mass density profile for the entire cluster, by 
combining the Subaru weak-lensing measurements with the 
inner strong-lensing information from deep, high-resolution 
HST/ACS/NIC3 observations (Zitrin et al. 2009b). 

5.1. One Dimensional Analysis 
HST/ACS/NIC3 Strong Lensing Constraints 

As demonstrated in Broadhurst et al. (2005b) and Umetsu 
& Broadhurst (2008) it is crucial to have information on the 
central mass distribution in order to derive useful constraints 
on the degree of concentration in the cluster mass distribution. 

To do this, we constrain the two NFW parameters from 
fitting to the combined weak and strong lensing data: 

X^ = xk + XsL> (32) 
where the Xwl weak lensing is defined by 

xk = E(<'^')-^0 . (33) 

with ki = k(0i) being the NFW model prediction for the lens- 
ing convergence at 0i and (C~'),j being the inverse of the 
bin-to-bin error covariance matrix for the one-dimensional 
mass reconstruction. The Subaru outer k profile is given 
in 10 logarithmically-spaced bins in radius 9 = [0.66', 10']. 
For the strong-lensing data, we utilize the improved strong- 
lens model of Zitrin et al. (2009b), well constrained by 33 
multiply-lensed images. With this model, we calculate the 
iimer k profile around the dark-matter center (§4.2) in 16 
Unearly-spaced radial bins spanning from ~ 0.13' to ~ 0.80', 
and the amplitude is scaled to the mean depth ~ 0.612 
of our blue-nred background sample (Table 3). For a joint 
fit, we exclude the strong-lensing data points at radii over- 
lapping with the Subaru data, yielding 12 independent data 
points as strong-lensing constraints for out joint fit. Finally, 
the function for the strong-lensing constraints is expressed 
as 

XL = E ^2 (34) 

where k, = k(0i) is the model prediction of the NFW halo for 
the ;th bin, and ct, is the Icr error for (k,) = {K(0i)); the bin 
width of the inner k profile is sufficiently broad to ensure that 
the errors between different bins are independent. By com- 
bining the full lensing constraints from the ACS/NIC3 and 
Subaru observations (22 data points), we can trace the cluster 
mass distribution over a large range in amplitude k ~ [10"^, 1] 
and in projected radius R = Dd9c^ [40, 2300] kpc h~^. 



In Figure 14 we show, for the entire cluster, the radial 
profile of the dimensionless surface mass density k{0) = 
E„(0)/Scrit from our combined strong- and weak-lensing ob- 
servations, where all of the radial profiles are scaled to a fidu- 
cial redshift of Zj = 1. For comparison purposes, the inner 
K profile {triangles) is shown to the maximum radius ~ 0.83' 
probed by the 33 multiply-lensed images, allowing for a direct 
comparison of the strong and weak lensing measurements in 
the overlapping region. This comparison convincingly shows 
that our strong and weak lensing measurements are fully con- 
sistent in the overlapping region. 

The resulting constraints on the NFW model parameters 
and the predicted Einstein radius 0-e (at Zs = 1 .675) are shown 
in Table 4. We show in Figure 15 the 68.3%, 95.4%, and 
99.7% confidence levels in the Cyir-Myir plane, estimated from 
Ax^ = 2.3,6.17, and 11.8, respectively, for each of the Subaru 
{dashed contours), ACS/N1C3 {solid contours), and joint Sub- 
aru and ACS/N1C3 {filled gray areas) k data sets. Apparently 
the constraints are strongly degenerate when only the inner 
or outer k profile is included in the fits. The virial mass Myu- 
is well constrained by the Subaru data alone, while the Sub- 
aru constraint on Cyir is rather weak, allowing a wide range 
of the concentration parameter, Cvir (see §4.4). On the other 
hand, the inner k profile from ACS/NIC3 observations probes 
up to ^ ~ 0.8', or the projected radius of /? ~ 0.18Mpc/i"' at 
the cluster redshift, which however is only about one-tenth of 
the cluster virial radius inferred from our full lensing analysis 
of Subaru and ACS/N1C3 data, r^i, ~ 1.8Mpc/z"', resulting 
in a rather weak constraint on Myir- Combining complemen- 
tary strong- and weak-lensing information significantly nar- 
rows down the statistical uncertainties on the NFW parame- 
ters, placing stringent constraints on the entire mass profile: 
= 1.15!!J:|5« X IO'^Mq/!-' and c^r = 9.21\i (xL/dof = 
4.1/20). All the sets of NFW models considered here are 
consistent with each other within the statistical uncertainty, 
and properly reproduce the observed location of the Einstein 
radius (see Table 4). 

Our high-quality lensing data, covering the entire cluster, 
allow us to place useful constraints on the gNFW structure pa- 
rameters, namely the central cusp slope a as well as the NFW 
virial mass and concentration parameters. Using our full 
lensing constraints, we obtain the best-fit gNFW model with 
the following parameters (see also Table 4): a = 0.04'!^^'^, 
Mv,r = 1 .06!;j;|^ X IO'^Mq r\ and c-2 = 9.8 ± 1 .2. The result- 
ing best-fit K profile from our full lensing analysis is shown in 
Figure 14 as a thin solid {gray) curve, along with the best- fit 
NFW profile {a = 1). Our combined weak and strong lensing 
data of CL0024H-1654 favor a shallower inner density slope 
with a < 0.97 (68.3% CL). The two-dimensional marginal- 
ized constraints (68.3%, 95.4%, and 99.7% CL) on {M^t,a) 
and {c-2,a) are shown in Figure 16. We note that the devi- 
ation in the inner density profile between the best-fit gNFW 
and NFW models becomes significant only below the inner- 
most radius of our lensing data, corresponding to the size of 
the radial critical curve. This results in a rather poor constraint 
on the central cusp slope (cf. Newman et al. 2009). 

Einstein-Radius Constraint 

As a model-independent constraint, we can utilize the ob- 
served location of tangential critical curves traced by gi- 
ant arcs and multiply-lensed images of background galaxies, 
defining an approximate Einstein radius, 6e (Broadhurst et al. 
2005a; Oguri & Blandford 2009; Zitrin et al. 2009b; Richard 
et al. 2009). For an axiaUy- symmetric lens, the Einstein- 
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radius constraint is written as k(< 6e) = 1, or ^+(6'e) = 1 (see 
equation [25]), corresponding to the maximum distortion, and 
provides an integrated constraint on the inner mass profile 
interior to 6e. More generally, an effective Einstein radius 
can be defined by axially averaging the projected surface- 
mass density, which itself is well determined when a large 
number of constraints are available (Broadhurst & Barkana 
2008; Richard et al. 2009). The Einstein-radius constraint for 
C10024H-1654 is shown in Figure 11 as the innermost data 
point (see also Figure 12). Here we follow the method de- 
scribed in Umetsu & Broadhurst (2008, §5.4.2) for incorpo- 
rating the inner Einstein-radius information into lensing con- 
straints: See also Oguri et al. (2009). 

We constrain the NEW model parameters (cvir,Mvir) by 
combining the Subaru g+ profile and the Einstein-radius in- 
formation. We define the function for combined lens- 
distortion and Einstein-radius constraints by^^ 



X 



(35) 



where the first term is the x^-function for the Subaru tan- 
gential shear measurements (g+^j) = {g+iOi)) and the second 
term is the function for the Einstein radius constraint; 
g+j = g+ifii) is the NEW model prediction for the reduced tan- 
gential shear bX9 = 6i calculated for the blue+red background 
sample (see §3.3), g+£ = g+(OE,ZE) is the model prediction 
for the reduced tangential shear at 6* = 6e, evaluated at the 
arc redshift, Zs = Ze- Following Zitrin et al. (2009b), we take 
9e = 30" with an rms error of ag^ = 3", corresponding to the 
observed 5-image system at ze = 1-675. We then propagate 
this error to g+ as ct+^e = <^0E\dg+/d6E\e=eE — 0-2 (see Umetsu 
& Broadhurst 2008). Similarly, one can combine a k profile 
reconstructed from lens distortion measurements with irmer 
Einstein-radius information (see §5.4.2 of Umetsu & Broad- 
hurst 2008). 

We show in Figure 17 the resulting constraints in the Cyir- 
Myir plane obtained from the Subaru g+ profile with (filled 
gray areas) and without (dashed contours) the inner Einstein- 
radius information, together with the NEW Cyir-Mvi,- rela- 
tion for the observed Einstein-radius constraint (9e = 30" at 
Ze = 1.675). The constraints from strong and weak lens- 
ing are fairly consistent with each other, showing similar 
degeneracy directions in the Cyir-Myir plane. Our joint fit 
to the Subaru g+ profile and the inner Einstein-radius con- 
straint tightly constrains the NEW model parameters, Mvir = 
1.19tix IO'^Mq/z-i andcvir = 8.6!j j (xL/dof =4.9/9), in 
good agreement with those from a joint fit to the ACS/NIC3 
and Subaru mass profiles (§5.1.0). 

5.2. Two-Dimensional Analysis: Two-Component Lens 
Model 

Our one-dimensional treatment thus far does not take into 
account the effect of the NW clump located at a projected dis- 
tance of ~ 3' (§4.1). Here we aim to account for this effect 
by two-dimensional lens modeUng with a multi-component 
shear fitting technique (Kneib et al. 2003), in conjunction with 
the inner strong-lensing constraints on hi(6). This method uti- 
lizes unbinned distortion measurements for individual galax- 

Strictly speaking, the rms dispersion cr, for the distortion measurement 
Si = g(Oi) is given as tr, ,(1 -\g(Oj)\^) in the subcritical, nonUnear regime 
(Schneider et al. 2000). We however neglect this nonlinear correction for the 
shear dispersion, and adopt a weighting scheme as described in §3.2.0. 



ies, and hence does not require any binning of distortion data 
(cf. §5.1). 

Assuming particular mass profiles for lens components, the 
model reduced-shear field g(9) = g\(6) + ig2(6) is simply cal- 
culated as 



(36) 



where Ni is the number of lens halo components. Kneib et al. 
(2003) analyzed a sparse- sampled mosaic of 2-band (F450W, 
F814W) WFPC2 observations in the C10024+1654 field, and 
found that two lens-components are necessary to match their 
WFPC2 lens-distortion data (i.e., A^; = 2), responsible for the 
central and NW clumps in projection space. 

Following Kneib et al. (2003), we use a circularly- 
symmetric NEW profile to describe the projected lensing 
fields of the central component, which has been shown to be 
a good approximation from our one-dimensional full lensing 
analysis (§5.1). Further, we assume that the central compo- 
nent is responsible for the central strong-lensing constraints 
on K(d) derived from the ACS/NIC3 observations. For the 
NW component, we use a truncated form of the NEW profile 
(tNEW, see Takada & Jain 2003), which approximates a struc- 
ture of a stripped halo by a sharp truncation at the halo virial 
radius. The location of the central component is fixed at the 
dark-matter center of mass (§4.2), around which the inner k 
profile (at 0.17' <e< 0.80' with 16 bins) is defined (§5.1.0). 
Further, the location of the NW component is fixed at the ob- 
served density peak position of the NW clump of green cluster 
galaxies (§4.1). 

We constrain two sets of NEW model parameters (Myir, Cyir) 
for the central and NW lens components by minimizing the 
total x^ function for our combined two-dimensional weak- 
lensing distortion and central strong-lensing constraints. The 
total x^ function is given as equation (32), but with the fol- 
lowing x^ function for Subaru weak lensing: 



2 _sr^ sr^ igaiOd-gaiOi)] 
XWL,2D - 2^ 



i=l a=l,2 



(37) 



where the index ; runs over all objects in our blue+red sam- 
ple, but excluding those at radii overlapping with the in- 
ner strong-lensing data (^^ ^ 0.80'), and (Jg,a(6i) (a =1,2) 
is the rms error for the real/imaginary component of the 
ith reduced-shear measurement gai^i), which we take as 
i7g,i(0i) = Og^2(0i) = ag(ei)/V2. We have 11647 such usable 
objects in the blue+red sample, i.e., a total of 2 x 11647 in- 
dependent measurements for the spin-2 distortion field. Thus, 
we have a total of 23310 joint constraints from strong and 
weak lensing, and 23306 dof. 

Table 5 lists the resulting best- fitting parameters of our two- 
component lens mass model. The two-component lens mass 
model provides an acceptable fit with the minimized x^ value 
of Xmin/dof = 1.577, and with the best-fit NEW parameters 
for the central component, Myk = (1 . 1 1 ± 0. 18) x 10'%© 
and Cvir = 8.1 ± 1.2, fairly consistent with those from the cor- 
responding one-dimensional analysis (§5.1). For the NW 
halo component, we find a best-fit set of NEW parameters, 
Mv„ = (1.28±0.51) X IO'^Mq/i-' and Cvir = 5.0±3.5, with 
the virial radius (and hence the truncation radius), ryir = 
(0.83 ±0.1 l)Mpc/j~\ corresponding to the angular radius of 
6lvir = 3.7' ± 0.5'. It is found that the best-fitting NEW pa- 
rameters obtained with and without the truncation at the virial 
radius agree to within 2% for the NW clump. 
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Our results can be directly compared with those from Kneib 
et al. (2003), who obtained M200 = 1.7![J| x lO^^^Qr^ and 
C200 = for the NW clump, where the quantities here are 
determined at r2oo corresponding to a mean interior over- 
density of 200, relative to the critical density Pcritfe) of the 
universe at the cluster redshift Zd = 0.395. These parame- 
ters can be translated into the corresponding virial param- 
eters, Mvir = 2.0:I^ Q X IO'^^Mq/i"' and Cyk = 5% consistent 
with our results within the statistical uncertainty. For the cen- 
tral component, on the other hand, Kneib et al. (2003) ob- 
tained an NFW model with very high concentration, M200 = 
A.G^j X 1O''*M0 /;"' and C200 = 22!^, or in terms of the virial 
parameters, Mvir = 4.3!^| x IQI'^Mq/j-' and Cyir = 26!^^ This 
high concentration may be due to the inclusion of the Einstein 
radius constraint in their fit to their outer weak lensing data 
which would otherwise be underestimated by weak lensing 
data alone (Umetsu & Broadhurst 2008, §5.5.2). It could also 
be partly explained by the irregular, non-axisymmetric mass 
distribution in the cluster core (Comerford et al. 2006; Comer- 
ford & Natarajan 2007). Simply adding the two virial masses 
of Kneib et al. (2003) yields Mot = 6.3!|:| x IO'^Mq r\ rep- 
resenting a large discrepancy with respect to our results. 

The radial mass profile of C10024+1654 has also been ex- 
amined by Hoekstra (2007) using ground-based weak-lensing 
shape measurements from CFHT/CFH12K R-bwd data. By 
fitting an NFW profile to the tangential distortion signal at 7? = 
[0.25, 1.5]Mpc/!~', Hoekstra (2007) obtained a virial mass of 
Mvir = 1 .76!^ ;g X 10'^Mq for this cluster, which is higher 
than, but consistent with, our results. 

5.3. Mass-to-Light Ratio 

Having obtained the radial mass density profile, we now 
turn to examine the cluster mass-to-Ught ratio (M/L) in 
a model independent approach. For this we utilize the 
weak lensing dilution method developed by Medezinski et al. 
(2007) (see also Medezinski et al. 2009), and derive the clus- 
ter luminosity density profile T^iiO) to large radius, with the 
advantage that no subtraction of far-field background counts 
is required. We weight each "green" galaxy flux Fi by its tan- 
gential distortion, g+j, and subtract this "^+-weighted" lumi- 
nosity contribution of each galaxy, which when averaged over 
the distribution will have zero contribution from the unlensed 
cluster members. This will account for any difference in the 
brightness distributions of the cluster members to that of the 
background galaxies, in particular the skewness of the cluster 
sequence to brighter magnitudes. The total flux of the cluster 
in the nth radial bin is then given as 

iebinn \S+,nl iGbinn 

where (g+^n) = (^^'(^n)) is the true background level of 
the tangential distortion, averaged over the nth radial bin, 
and (/3*^') and (Z?*^^') are the source-averaged distance ra- 
tios (see equation [7] and Table 3) for green and reference 
background samples, respectively (see Appendices A and 
B of Medezinski et al. (2007) for a derivation of this equa- 
tion). Here we take the blue-nred sample as our reference 
background. The flux is then translated to luminosity. First 
we calculate the absolute magnitude, 

M = m-5\ogddz)-K{z) + 5, (39) 

In general, background samples may contain foreground field galaxies 
with Zs < Zd and = 0. 



where d^iz) is the luminosity distance to the cluster, and the 
/T-correction, K{z), is evaluated for each radial bin according 
to its colour. We use the /?c-baiid data to calculate the cluster 
Ught profile. 

The results of the ^T-corrected cluster luminosity density 
measurements are shown in Figure 18, along with our lensing 
constraints on the surface mass density profile, Sm(^) {shaded 
gray area), converted into a luminosity density profile as- 
suming a constant M/Lr of 22>Qh{M /Lr)q, corresponding to 
the mean cluster mass-to-light ratio interior to 9 = 3' . There 
are two sets of S; profiles shown in Figure 18, namely with 
(squares) and without (circles) the contribution from the NW 
galaxy clump located at a projected distance of 6* ~ 3' (§4.1). 
The observed light profile, including the NW clump, closely 
resembles the mass profile derived from our joint weak and 
strong lensing analysis, and shows a shoulder feature at ^ ~ 
3'. This feature disappears in the NW-clump corrected light 
profile, and hence is caused by the excess luminosity due to 
the NW galaxy clump. The total luminosity of the NW clump 
is estimated as L(< V) = (3.84±0.20) x IO'^Lq/j-^ andL(< 
1.5') = (5.47 ±0.45) X IO^Lq/j-^ for aperfiires of radius 1' 
and 1 .5', respectively. Assuming the same mean mass-to-ratio 
of {M/Lr) = 230/!(M/Lr)q for the NW galaxy clump, we 
have projected mass estimates of M2d(< 1') = (0.86 ± 0.05) x 
lO'^M^ri and M2d« 1.5') = (1.23 ±0.10) x lO^^^or^; 
these values are consistent, within the errors, with the predic- 
tions of our NFW model constrained by the combined weak 
and strong lensing data. 

We show in Figure 19 our model-independent radial pro- 
file of the differential mass-to-hght ratio, SM(9)/6Lr(9) = 
T,m(9)/'Ei(6), obtained by dividing the non-parametric 
profile from our Subaru distortion data by the /?c-band S; pro- 
file. Our direct, model-independent approach yields larger er- 
rors than those from conventional parametric procedures as- 
suming a particular form (e.g., an NFW profile) for the un- 
derlying mass density profile. The resulting 5M/5Lr pro- 
file is noisy and consistent with a constant mass-to-light 
ratio, but shows a local peak around R ~ IMpc/z"' with 
M/Lr = 400!^^[;/z(M/Lr)q, and then falls off to larger ra- 
dius. A similar behavior has been convincingly shown by a 
recent lensing-based analysis (Medezinski et al. 2009) of sev- 
eral high-mass clusters (A1689, A1703, A370, RXJ1347-11) 
with Mvir^ IO'^Mq/j-i using deep multi-color Subaru im- 
ages. For these relaxed, high-mass clusters, it has been 
shown that the S; profiles all dechne smoothly as ~ in 
projection, whereas the profiles are well described by 
the continuously varying NFW profile, so that the SM/SL 
profile peaks around 20% of the virial radius in the range 
(3OO-5OO)M0/L0, and then steadily falls to a mean field 
value of ~ IOOM0/L0, consistent with careful dynamical 
work by Rines et al. (2000). In contrast, the M/Lr peak 
for C10024H-1654 is found to be at ^ 0.6rvir, much larger 
than those found for other massive, relaxed clusters with cen- 
trally peaked M/L profiles. This different radial trend of 
C10024H- 1654 may reflect a different dynamical state from re- 
laxed systems. 

6. MASS PROFILE AND COMPARISON 

6.1. Comparisons with Other Lensing Studies 

First we compare our results with those from previous lens- 
ing studies of C10024H-1654 to assess the consistency over 

Medezinski et al. (2009) used their best-fit NFW S„ profile to derive 
cluster mass-to-light ratio profiles. 
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different regimes of gravitational lensing (weak and strong 
regimes) and to highlight any potential problems or discrepan- 
cies arising from systematic effects, such as weak lensing di- 
lution due to contamination by unlensed foreground and clus- 
ter member galaxies, shear calibration biases, and spurious 
boundary effects. 

Tangential Distortion and Lensing Convergence Profiles 

In Figure 20 we compare our Subaru tangential-distortion 
profile {squares) with the results from different authors and 
different observations/methods: WFPC2 weak-lensing analy- 
sis of Kneib et al. (2003), ACS weak-lensing analysis of Jee 
et al. (2007), and ACS/NIC3 strong-lensing analysis of Zitrin 
et al. (2009b). For this comparison, we have converted the 
azimuthally-averaged convergence profile n{9) of Zitrin et al. 
(2009b) into a distortion profile {shaded region) for a fiducial 
source at = 1.3, roughly matching the mean depth of our 
Subaru blue-nred background sample (Table 3). 

Our Subaru and ACS/NIC3 results are in good agreement in 
the overlapping region (0.6' ^9^1'). Furthermore, a sim- 
ple extrapolation of the best-fitting NFW profile for the in- 
ner ACS/N1C3 observations {solid; see §5.1.0) fits well with 
the outer Subaru distortion information over a wide range 
of radius, but somewhat overpredicts the distortion profile at 
3' ^ ^ ^ 5', as our Subaru data prefer a slightly steeper pro- 
file. The triangles show the g+ profile as obtained from the 
ACS weak-lensing analysis of Jee et al. (2007), where the data 
are limited to the positive-parity region {0 ^ 40") in this com- 
parison. Their overall g+ profile is steeper than our combined 
ACS/N1C3 and Subaru results, but is roughly consistent out to 
^ ~ 1 .5', beyond which it drops more rapidly than our Subaru 
profile. The open circles show the g+ profile at 0.9' ^0 ^ 11' 
from the weak-lensing analysis by Kneib et al. (2003). The 
mean depth of their background sample iszj=1.15±0.3 from 
photometric redshift measurements in the Hubble Deep Field- 
South (Kneib et al. 2003). The flattened slope of Kneib et al. 
profile at 9^2' is likely due to contamination of the weak 
lensing signal by unlensed cluster member galaxies. Indeed, 
Kneib et al. (2003) made an empirical correction for this di- 
lution effect in their two-dimensional lens modeling. Never- 
theless, taking into account the difference in redshift depths, 
our Subaru and Kneib et al. results are in agreement at ^ 3', 
where the weak-lensing approximation is valid {\g\ ^0.1) and 
the dilution effect due to contamination by unlensed cluster 
galaxies is less significant. At the sensitivity and resolution of 
ground-based weak-lensing measurements, our Subaru distor- 
tion profile shows no evidence for the dip at ~ 1 .25', corre- 
sponding to the ringlike structure revealed by the ACS lensing 
analysis of Jee et al. (2007). 

Similarly, we compare in Figure 21 the radial profiles of 
lensing convergence k{9) = S„(^)/Scrit, or the dimensionless 
surface mass density, derived from the lensing studies shown 
in Figure 20. Here all of the k profiles are scaled to the same 
fiducial source redshift of Zs = 1 for direct comparison. As we 
have already shown in Figure 14, our weak and strong lens- 
ing results are fully consistent with each other in the overlap- 
ping region, and the combined k profile of the full lensing 
constraints, for the entire cluster, is well fitted with a sin- 
gle NFW profile (see §5.1.0). At larger radii of 9 ^ 2', our 
Subaru-derived k profile agrees well with the two-component 
NFW lens model of Kneib et al. (2003), constrained from their 
weak-lensing measurements combined with the inner strong- 
lensing information on the Einstein radius, whereas their k 
values at smaU radii, 0.3' ^9^2', are systematically lower 



than our Subaru results. On the other hand, the n profile re- 
constructed from the ACS lensing analysis of Jee et al. (2007) 
is in agreement with our Subaru and ACS/NIC3 results at 
radii 0.3' ^6 ^ 0.8', beyond which, however, the surface- 
mass density of Jee et al. (2007) stays almost constant at 
K ~ 0.5 out to the maximum radius of the ACS observations 
(^max ^ 370kpc/!"' at z = 0.395), and largely disagrees with 
both our and Kneib et al.'s results. The overall shallow pro- 
file obtained by Jee et al. (2007) could be due to the mass- 
sheet degeneracy (§2), as demonstrated by Liesenborgs et al. 
(2008b). 

Cumulative Projected Mass Profiles 

Now we derive a non-parametric projected mass profile 
M2d(< 9) = Tr{Dd9)^T,cntK{< 9) of C10024H-1654 from our 
combined k profile of the ACS/NIC3 and Subaru lensing ob- 
servations described in §5.1.0; there are a total of 12 and 10 
data points from the ACS/N1C3 and Subaru observations, re- 
spectively, over the radial range R ~ [40,23001kpc/!"'. We 
interpolate k{9) between measured points with a cubic spline 
function. We use Monte-Carlo techniques to estimate the 
confidence limits on M2d(< 9), by properly propagating the 
measurement errors encoded in the joint covariance matrix 
Cij = {SkiSkj), containing both diagonal and off-diagonal ele- 
ments, constructed from the combined ACS/N1C3 and Subaru 
K data. To do this, we Cholesky-decompose this covariance 
matrix as C = LL^, where i is a lower triangular matrix, and 
assign a random noise fiuctuation 5k{9i) in each radial bin by 
Ski = LijS^j with being drawn from the Gaussian distribu- 
tion with zero mean and unit variance (e.g., see Park et al. 
2003). We then generate 1000 random reaUzations of the k 
profile, and perform an identical analysis on the Monte-Carlo 
data sets. 

Figure 22 shows the resulting projected mass profile Mxd 
of C10024H-1654, along with the results from Kneib et al. 
(2003) and Jee et al. (2007). The resulting 68.3% confi- 
dence limits from our combined ACS/NIC3 and Subaru ob- 
servations are shown as the gray-shaded area, in good agree- 
ment with the results from Jee et al. (2007, Figure 12) at 
radii up to ^ I' {R ^ 220 kpc /;"'). Beyond this radius, how- 
ever, their mass profile shows a much steeper increase out 
to the boundaries of the ACS observations (see also Figure 
21), owing to the flat behavior of their k profile at 6* ^ 0.8' 
(i.e., M2D 9^), and thus disagrees with our mass profile. 
Beyond 9^ — 0.5' for the 5-image system (z, = 1.675), the 
profile from the best-fit two-component NFW model of 
Kneib et al. (2003) is systematically lower than our and Jee et 
al.'s mass profiles. Finally, all of these results are in close 
agreement at radii around the Einstein radius 9e — 30" at 
Zs = 1.675, tightly constrained by the 5-image system. From 
our non-parametric M2D profile, we find a projected mass of 
M2d(< 30") = (1.32±0.12) X lOi'^r', in excellent agreement 
with M2d(< 30") = (1.27 ±0.09) x IO^'^MqH-^ obtained by 
Jee et al. (2007). 

6.2. Mass Discrepancy 

A comparison of cluster mass estimates from X-ray, dy- 
namical, and strong/weak lensing observations is of particular 
importance to assess the vaUdity of fundamental assumptions 
made in various methods of mass determination. As com- 
pared with gravitational lensing, which is purely geometrical 
and sensitive to the projected mass along the line of sight, X- 
ray and dynamical methods infer the three-dimensional struc- 
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ture of clusters under the assumptions of symmetry and dy- 
namical equilibrium. In particular, the X-ray method can be 
affected strongly during mergers where the assumption of hy- 
drostatic equilibrium is no longer valid (Okabe & Umetsu 
2008). The redshift-space caustic method of Diaferio (1999), 
in contrast to the traditional Jeans approach, does not rely on 
the assumption of dynamical equilibrium, and can be used to 
estimate the cluster mass even in non-equilibrium regions. 

A detailed dynamical analysis by Czoske et al. (2001, 
2002) revealed a fairly complex, bimodal velocity structure in 
C10024-I-1654, with two distinct redshift components peaked 
at 2 = 0.395 (component A) and z = 0.381 (component B). The 
former is the primary cluster, and the latter is a foreground, 
loose aggregate of galaxies physically separated from the pri- 
mary component. They showed that the redshift distribution 
of cluster members in the central region is strongly skewed 
towards negative velocities. Czoske et al. (2002) suggested 
that the diffuse foreground component could be the result of 
a high-speed, head-on encounter of two similar-mass clus- 
ters with a merger axis very nearly parallel to the line of 
sight. In this context, the relative velocity of the two com- 
ponents implied by the redshift difference is v ^ 3000km s~'. 
With A'-body simulations of two colliding CDM halos with 
a mass ratio of about 2:1, they also demonstrated that such 
a scenario can explain the observed peculiar redshift distri- 
bution in C10024H-1654. In their simulations they adopted 
Mam, = 9.5 x IO^^Mq/;-' and Mbm, = 5.0 x IO^^Mq/z"' as 
initial virial masses for the components A and B, respectively. 
It is interesting to note that the sum of the two virial masses, 
Mvir = 1 .45 X 1O'^M0 , is close to our la upper Umit on the 
total virial mass Mvir — 1.4 x lO'^M©/;"', and is consistent 
with the total projected mass estimate M2d(< ''vir) = 1 -SS^y ^j 
for the entire cluster system, including any currently unbound 
material beyond the virial radius. 

The collision scenario of Czoske et al. (2002) has been also 
supported by the presence of a ringlike structure in the pro- 
jected dark matter distribution revealed by the ACS lensing 
analysis of Jee et al. (2007). They speculated from the dis- 
sipation of the shock in the Chandra and XMM-Newton X- 
ray observations and the size of their observed ringlike struc- 
ture (R ~ 280kpc/j~^) in their mass map that the two clusters 
in C10024H-1654 would have undergone the first core passage 
about 1-2 Gyr ago, and that the dense gas cores of the two 
clusters would have survived the collision with a distinct sep- 
aration. If the cluster is not considered to represent a single 
relaxed system, but is a superposition of two separated clus- 
ters, then the single-component assumption should lead to a 
substantial underestimation of the X-ray and dynamical mass, 
whereas the lensing mass should be the sum of the two cluster 
components. 

To highlight this problem, we now compare our results with 

integrated three-dimensional mass profiles Mj,o(< r) derived 
from X-ray and dynamical, as well as lensing, observations 
with the single cluster assumption. Here we deproject the 
two-dimensional mass profile and obtain a non-parametric 
M3D profile simply assuming spherical symmetry, following 
the method introduced by Broadhurst & Barkana (2008). This 
method is based on the fact that the surface-mass density 
E„,(/?) is related to the three-dimensional mass density p(r) 
by an Abel integral transform; or equivalently, one finds that 
the three-dimensional mass M3d(< r) out to spherical radius 



r is written in terms of T.,„{R) as 

M3D« r) = 2ttJ^ dRRE„iR)-4j^ dRRf (^^^ S„(^), 

" ' (40) 

where f(x) = (x^ - ly^'^ -taiC^x^ - l)'^'^ (Broadhurst & 
Barkana 2008). To perform the integrals in equation (40), 
we assume smoothness of Sm(^) and use the interpolation 
method described in §6.1.0. We then extrapolate inward as- 
suming I^miR) = const, from the innermost point and outward 
assuming S;„(/?) oc R~^ from the outermost point, where these 
exponents are motivated by the NFW profile, which provides 
an excellent description of our data (§5.1.0). We find that 
even varying these exponents by ± 1 would change the mass 
profile at all relevant radii by only ^ 0.3%, which is negli- 
gible compared to the effect of the measurement errors. Fi- 
nally we propagate errors on T,m{R) using the same Monte- 
Carlo method as in §6.1.0. This deprojection method allows 
us to derive in a non-parametric way three-dimensional virial 
quantities (rvir,Mvir) and values of Ma = M3d(< ta) within a 
sphere of a fixed mean interior overdensity A with respect to 
the critical density of the universe at the cluster redshift. Table 
6 gives a summary of the cluster mass estimates M3d(< ta) 
from our non-parametric deprojection analysis. 

In Figure 23, the resulting deprojected mass profile M3d(< 
r) and 68.3% confidence level region are shown as the solid 
line and gray-shaded area, respectively, and compared with 
results from other studies. Also indicated in Figure 23 is 
the location of the virial radius ryir = 1.75^q jj Mpc/i"', de- 
rived from the non-parametric deprojected mass profile, cor- 
responding to the virial mass ofMvir= 1.21^q22 x IO'^Mq/z"', 
which is in excellent agreement with those from our best-fit 
NFW models, M™ = (0.9-1.3) x IO^^Mq h-^ (la). Our non- 
parametric estimate of M500 is M500 = 7.2^2 ^ IO'^^Mq/i"' 
within the spherical radius r^oo = 0.93^g jqMpc/z"', in good 
agreement with the model-independent mass estimate M500 = 
(7.4 ± 2.4) X IO'-^Mq r 1 by Hoekstra (2007). 

Clearly, the comparisons in Figure 23 exhibit a large dis- 
crepancy between our lensing and the X-ray/dynamical re- 
sults based on the single cluster assumption, indicating a sub- 
stantial contribution from additional mass components. The 
dashed curve shown is the spherical NFW mass profile from 
the best-fit lens model of Kneib et al. (2003), where only the 
central halo component has been included. The X-ray derived 
mass profile (dotted curve) of Ota et al. (2004) is taken from 
Figure 1 of Diaferio et al. (2005). Overall, the caustic mass 
estimates (triangles) of Diaferio et al. (2005) he between the 
results of Kneib et al. (2003) and Ota et al. (2004). 

OurM3D profile is only marginally consistent with the caus- 
tic mass estimates of Diaferio et al. (2005) at r ^ 0.7Mpc/!"\ 
beyond which the mass discrepancy increases with increasing 
radius. The outermost radius probed by Diaferio et al. (2005) 
is close to our estimated virial radius (Table 6), at which the 
lensing to caustic mass ratio is about 3.0 ± 1.5. This mass 
discrepancy could be largely reconciled within the collision 
scenario of Czoske et al. (2002), in which C10024H-1654 is 
the result of a high-speed collision of two massive clusters 
along the line of sight. Simply assuming that the total virial 
mass of the whole system is conserved before and after the 
collision, with a fiducial mass ratio of 2: 1 (Czoske et al. 2002; 
Jee et al. 2007; Zu Hone et al. 2009b,a), the total mass of the 

" This integral transformation lias an integrable singularity at the lower 
limit of the second integrand (R = r), which can be readily avoided by a suit- 
able coordinate transformation. 
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primacy cluster (A), excluding the NW substructure, is then 

marginally consistent with 
the caustic mass estimate within the uncertainty. Further, such 
a collision may disrupt the dark matter from the cluster cores, 
and eject a substantial amount of mass into cluster outskirts, 
while the velocity dispersion at large radii can be significantly 
reduced (by 10%- 15%, Czoske et al. 2002); a detailed dis- 
cussion will be presented in §7.1. This may lead to a system- 
atic mass discrepancy between lensing and caustic methods at 
cluster outskirts. 

Under the hypothesis of hydrostatic equilibrium, Ota 
et al. (2004) derived a spherical hydrostatic mass of M200 = 
3.2+',;^ X lO"*M0/z-' within a sphere of raoo ^ LOMpc/z"' 
(^^200 = 4.3')- At this radius, our derpojected mass estimate 
is M3D = 7.4;^;j 4 X 1O'^M0 /;"'. This comparison yields a lens- 
ing to X-ray mass ratio of 2.3 ±0.7 at r ^ 1 Mpc/z"'. Inter- 
estingly, the Chandra and XMM-Newton observations of this 
cluster revealed that the X-ray surface brightness distribution 
is better fitted by a superposition of two ICG models rather 
than one (Ota et al. 2004; Zhang et al. 2005). Jee et al. (2007) 
demonstrated that assuming a superposition of two isothermal 
models for the X-ray brightness distribution can indeed bring 
the Chandra data into closer agreement with the lensing mass 
estimates, Ml^{< 30") ~ 1.1 x IQ^'^MqU-^ (see §6.1.0). 

7. DISCUSSION 
7.1. Post-Collision Scenario 

Our full lensing analysis of joint Subaru and ACS/NIC3 
observations has shown that the radial structure of the pro- 
jected mass in C10024H-1654 is consistent with a continuously 
steepening density profile, with only a minor contribution 
from the NW substructure, and well described by the gen- 
eral NFW profile, as previously found for a number of re- 
laxed, massive clusters (Broadhurst et al. 2005b; Umetsu & 
Broadhurst 2008; Umetsu et al. 2009; Okabe et al. 2009), 
expected for collisionless CDM halos. Mass and light in 
the cluster are similarly distributed and exhibit a prominent 
peak at the cluster center in projection space. All these ob- 
served features might lead one to conclude that C10024H-1654 
is a fairly relaxed, massive cluster with a total virial mass of 
Mvir = (1 .2 ± 0.2) X lO'^M© /;"', but with noticeable substruc- 
ture of a fractional mass (1 1 ± 4)%, located at a projected ra- 
dius /? ~ 700 kpcr'. 

Several independent Unes of evidence, however, suggest 
that the cluster C10024H-1654 is the result of a high- velocity, 
near head-on collision between two similar-mass clusters 
occurring along the line of sight, viewed approximately 
2-3 Gyr after impact: 

(i) Dynamical data 

The first piece of evidence is the existence of a bimodal 
redshift distribution in the system (Czoske et al. 2002). In 
addition, the central velocity distribution of cluster members 
is strongly skewed negative with a long tail towards the 
foreground, secondary component, suggesting strong recent 
gravitational interaction (Czoske et al. 2002). 

(ii) X-ray emission features 

Next, anomalously low levels of X-ray emission and tem- 
perature (Ota et al. 2004; Zhang et al. 2005), relative to 
the standard X-ray observable and mass relations, may 



also point to interaction, and an incomplete merger. In a 
period of 1 - 3 Gyr ( ^ Tsc) after the first passage of a close 
encounter between two massive clusters, the shock-heated 
gas associated with this whole system could be in a much 
expanded phase with correspondingly lower X-ray luminosity 
and temperature (Ricker & Sarazin 2001). In this scenario, 
the dense gas cores of the two clusters are thought to have 
survived the collision with a distinct separation, and have not 
settled down to a single system (Jee et al. 2007). This is sup- 
ported by the high-resolution Chandra/XMM-Newton X-ray 
data which favor a superposition of two gas components over 
a single one (Ota et al. 2004; Jee et al. 2007). Recently, such 
X-ray brightness features observed in the cluster have been 
remarkably reproduced by a high-resolution simulation of a 
high- velocity head-on colUsion of 2: 1 mass-ratio clusters (Zu 
Hone et al. 2009b), motivated by the scenario of Czoske et al. 
(2002). In particular, it has been shown that a hydrodynamic 
mass estimate in a post-collision state, assuming one cluster 
component, underestimates the actual mass by a factor of 
2-3 (Zu Hone et al. 2009b), where the X-ray bias is found 
to be due mainly to the projection of denser, colder gas along 
the line of sight of hotter gas in the main cluster. 

(iii) Mass discrepancy between lensing and caustic methods 

Similarly, a mass discrepancy between lensing and caustic 
methods found in cluster outskirts (§6.2) may indicate the 
existence of an extended, diffuse mass component with 
lowered velocity dispersions due to evaporation of faster 
moving dark-matter particles and galaxies (Ricker & Sarazin 
2001; Czoske etal. 2002). 

The collision scenario of Czoske et al. (2002) appears to be 
consistent, at least qualitatively, with all these observations, 
and provides a feasible explanation for the long-standing puz- 
zle on the mass discrepancies and anomalies found in this 
cluster (§6.2). Nevertheless, a detailed and quantitative un- 
derstanding of the observed phenomena and the physical con- 
ditions of the hot gas and dark matter has yet to be obtained. 
Recent A^-body/hydrodynamic simulations of Zu Hone et al. 
(2009b) have been successful at reproducing qualitatively the 
observed features of X-ray brightness and temperature struc- 
ture in C10024H-1654. Their simulated X-ray observations 
show that the observed X-ray spectral temperatures can be 
lower than the actual temperature of the hotter gas by a factor 
of ~ 2 due to the superposition of denser, colder gas along the 
line of sight. However, the resulting X-ray temperature from 
their mock observations is only T ^ 2.6-2.8keV, which is ~ 
60% of the observed spectral temperature of Tx = 4.47^y 54 by 
Ota et al. (2004). Without resorting to further ad hoc explana- 
tion, this quantitative disagreement, in the context of the col- 
lision scenario of Czoske et al. (2002), indicates that the ini- 
tial conditions for their simulation are different from the pre- 
merger conditions in the real cluster (Zu Hone et al. 2009b). 
In particular, the sum of the initial masses of two clusters 
adopted by Zu Hone et al. (2009b) is M200 = 9 x 10"Mq 
{h = 0.5,f2m = 1 in their simulations) with a mass ratio of 
Ma 200 /Mb 200 = 2, and is ~ 60% of our lensing mass estimate 
(Table 6). 

Gravitational lensing, on the other hand, can provide es- 
sential information to constrain the total mass of the whole 
system in the pre-merger state. Our comprehensive lensing 
analysis provides a model-independent constraint of M2d(< 
rvir) = (1.36!^j^) X IQ^^MqU-^ for the projected mass of the 
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whole system, including not only the primary and secondary 
cluster components, but also any currently unbound material 
beyond the virial radius. If the collision is occurring along our 
line of sight, then the total projected mass associated with this 
system is essentially conserved before and after the colUsion. 
Thus, a constraint on M2d(< ''vir) can be regarded as a crude 
upper limit to the sum of the two pre-merger cluster masses 
when designing simulations to explore this system. 

The contrast between the largely monolithic lensing based 
mass profile and the bimodal post-merger interpretation of the 
dynamical and X-ray data may be largely reconciled when 
considering the radial dependence of the relaxation timescale. 
If we take the density at a radius of r = 1 Mpc where the mass 
density we estimate is p ^ 4 x 10"^^ gem"-', corresponds to 
a dynamical time of tdyn ~ 1 Gyr, shorter than or compara- 
ble to the estimated time since the merger occurred — which 
must be larger than ^1-2 Gyr so that the hot gas shock is 
no longer present and less than about 3 Gry so that the gas 
properties are still noticeably anomalous: This means that the 
central region, interior to 1 Mpc, can be expected to be well 
relaxed at this stage of the merger with a symmetric inner po- 
tential, which is born out on inspection of numerical simula- 
tions, in particular Ricker & Sarazin (2001). The outskirts, 
however, will not be relaxed for example at 2 Mpc where the 
mass density is p ~ 7 x 10~^^gcm~^, corresponding to a dy- 
namical time of tdyn ~ 4 Gyr, so that we can expect here to 
see large scale velocity bimodality remaining. However, the 
sum of two roughly aligned, unrelaxed diffuse outer halos of 
the two merging clusters will not be easily distinguished by 
lensing in projection. 

Recent results of a weak lensing analysis of Jee et al. (2007) 
reveal a ringlike structure in the central projected mass distri- 
bution. They interpreted this is the result of a high-speed col- 
lision between two dark-matter halos along the line of sight. 
They demonstrated that such ring features could be repro- 
duced in a simulation of a high-speed colUsion of two pure 
dark-matter halos. More recently, however, detailed A^-body 
simulations show that such features require significant fine- 
tuning of initial conditions, such that the initial particle ve- 
locity distribution is purely circular (Zu Hone et al. 2009a), 
which however is unUkely in the colUsionless CDM model. 
Our joint weak and strong lensing analysis shows no evidence 
for the presence of ringlike features at the sensitivity and reso- 
lution of ground-based weak lensing observations. Except for 
the presence of ringlike features, our joint distortion profile 
is in overall agreement with the results from a weak lensing 
analysis of Jee et al. (2007). On the other hand, our self- 
consistent joint K profile, with a continuously declining radial 
trend, is in disagreement with the results of Jee et al. (2007). 
These discrepancies between our results and those reported 
by Jee et al. (2007) may be explained by a combination of 
the mass-sheet degeneracy (§6.1.0) and the monopole degen- 
eracy, as suggested by Liesenborgs et al. (2008a,b). 

7.2. Implications for the Yx -Mtot Relation 

Recent cosmological hydrodynamicW-body simulations 

(Kravtsov et al. 2006) suggest that a simple product of the X- 
ray spectroscopic temperature and the X-ray derived hot-gas 
mass, namely, Yx = T^m x Mgas,5oo, is a robust X-ray indicator 
of the total cluster mass Mtot, being insensitive to the cluster 
dynamical state. The primary reason for the tight Ix-Mtot 
relation is that this X-ray mass proxy is directly related to the 
total thermal energy of the ICG, and thus to the depth of the 
cluster potential well. As a result, Yx is not strongly disturbed 



even by cluster mergers, unlike other X-ray mass proxies, Tx 
or X-ray luminosity (Ricker & Sarazin 2001). In addition, the 
stabihty of Yx can be explained by that fractional deviations 
of the average gas temperature Tsoo and X-ray derived Mgas,5oo 
tend to be anticorrelated (see Figure 5 of Kravtsov et al. 2006), 
and hence partially cancelled out in the product Yx. 

Here we use the Chandra X-ray data presented by Ota 
et al. (2004) to derive the value of Yx = Mgas.sooTsoo for 
C10024H-1654, assuming a single cluster component that rep- 
resents the primary cluster A. The outer radius for integration 
of the Chandra X-ray spectrum is only rout — O.Srsoo, where 
rsoo is based on the virial mass Mamt of the primary cluster 
component, MA,yii ~ 7.3 x IO^'^MqIi'^, obtained assuming a 
fiducial 2:1 mass-ratio collision scenario (§6.2). We estimate 
the value of Tsoo from the measured value Tx = T( ^ O.Srsoo) 
by using an empirical T^qo /Tx-Tx relation found by Vikhlinin 
et al. (2009, see equation |5] and Figure 6). The result- 
ing Yx from the Chandra X-ray observations is Fx (Chandra) = 
(3.8 ±2.3) x lOi^MokeV/ji/l This Chandra measurement of 
Yx is compared with that expected from the self-similar scal- 
ing relation Yx = ix(A^500,z) calibrated by detailed Chandra 
observations of a sample of nearby, relaxed clusters (Vikhlinin 
et al. 2009). Using this i^- A/500 relation, we obtain the 
value of Yx, predicted for the primary cluster component A, 
as Yx(Ma,5Qo) = (1.9 ±0.8) x lO'^MekeV/z'/^. Hence, the 
X-ray Yx measurement only accounts for ^ 35% (la) of 
the self-similar prediction based on the lensing mass esti- 
mate: Fx(Chandra)/Fx(MA,5oo) = 0.20 ±0.15. This discrep- 
ancy could be partially reconciled by the superposition effects 
of X-ray emission found by Zu Hone et al. (2009b). 

Another plausible mechanism to account for this apparent 
discrepancy is an adiabatic cooUng in an adiabatic expansion 
process. If the ICG is undergoing an adiabatic process, whose 
timescale is of the order of Tsc and much shorter than the cool- 
ing timescale over the cluster radius (e.g., rjoo), then the adia- 
batic condition, Tn^~^'^ = const, with 'jad being the adiabatic 
index of the ICG, relates the fluctuations in gas temperature 
and density as ST/T {2/3)Sn/n with jad = 5/3 (for non- 
relativistic, ideal gas). In X-ray observations, we measure the 
gas mass Mgas at a certain fixed radius (for example, provided 
by lensing), rsoo, so that Sn/n ^ <5Mgas,50o/A^gas,500- Since 
the adiabatic fluctuations 6T and 6n are positively correlated, 
the deviation in their product, Yx, is enhanced in an adiabatic 
phase: 

SYx _ (5Mgas,500 ^ STsQQ _ 5 6T500 ^^^^ 
Yx A/gas,500 7500 2 7500 

During the post-shock adiabatic expansion phase (t ^ Tsc), 
the fractional decrease in the average gas temperature could 
reach -ST /T ^ 20% at maximum for an head-on collision of 
two similar-mass clusters (Ricker & Sarazin 2001), suggest- 
ing that the maximum fractional decrease in the Yx is of the 
order of -SYx /Yx 50%, provided that the cluster is caught 
in such an adiabatic expansion phase. 

7.3. Numerical Simulations 

We evaluate this merger scenario based on numerical simu- 
lations of cluster collisions using FLASH3, an Eurelian three- 
dimensional hydrodynamics/A'-body simulation code (Fryx- 
eU et al. 2000). FLASH uses the Piecewise-ParaboUc Mefliod 
(PPM) of Colella & Woodward (1984) to solve the equations 
of hydrodynamics, and a particle-mesh method to solve for 
the gravitational forces between particles in the A^-body mod- 
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ule. The gravitational potential is calculated using a multi- 
grid solver (Ricker 2008). FLASH uses adaptive mesh refine- 
ment (AMR) with a tree-based data structure allowing recur- 
sive grid refinements on a cell-by-cell basis. With the AMR 
technique, our simulations achieve a maximum resolution of 
12.7 kpc. 

As pre-collision conditions, we assume two equilibrium 
NFW halos with a mass ratio of 2:1 following Czoske et al. 
(2002) and Jee et al. (2007). The sum of the virial masses 
of the two clusters is set to Myk = 1.1 x 1O'^M0 /z"' ~ 1 .57 x 
lO^^M© based on the results of our full lensing analysis, so 
that the initial virial masses for the primary (A) and sec- 
ondary (B) clusters are M^^vir = 1.05 x lO'^M© and MB,m = 
0.52 X 1O^^M0, with virial radii of rA,vir = 2.1Mpc and rB,m = 
1.7Mpc, respectively. The concentration parameters for the 
two clusters are set to c^^vir = 4.2 and c^ vir = 4.8, roughly 
matching the median Cvir(Mvir,z) relation for relaxed CDM 
halos in the WMAP five-year cosmology (Duffy et al. 2008). 
The clusters are initialized so that their centers are separated 
by the sum of their virial radii with an initial relative velocity 
of V = 3000kms-i (Czoske et al. 2002). We fix the gas mass 
fraction within rvir to /gas = 0.14, close to the values obtained 
from recent X-ray and multi- wavelength observations of mas- 
sive clusters (Vikhlinin et al. 2009; Umetsu et al. 2009). For 
simpHcity, we ignore the expansion of the universe during the 
collision. For the gas component of each cluster, we assume a 
(3 model in hydrostatic equilibrium with the gravitational po- 
tential of the cluster, with core radius = 0.08rvir and /3 = 1, 
motivated by results of our analysis of galaxy clusters drawn 
from high-resolution cosmological simulations (S. Molnar et 
al. 2010, submitted to ApJ). For the primary and secondary 
cluster, the maximum values of the three-dimensional tem- 
perature T^D are found to be 8.6 keV and 5.8 keV, respec- 
tively, in the radial range 0.15r5oo- ?"500. We use the local 
MaxweUian approximation for the local velocity distribution, 
and determine the dispersion from the Jeans equation assum- 
ing isotropic velocity distribution as a function of radius. The 
gas is simply assumed to have no bulk velocity relative to their 
respective NFW halos; the details of the simulations will be 
presented elsewhere (S. Molnar et al., in preparation). In the 
present study, we consider a head-on collision with zero im- 
pact parameter as done by Zu Hone et al. (2009b). 

We find that during the first core passage, which occurs at 
about f = 1 Gyr after the beginning of the colUsion, the in- 
faUing subcluster gas core pushes the main cluster gas core 
out of its equilibrium position. After the first core passage, 
the hot gas of the main cluster core, mixed with part of the 
gas stripped off the subcluster, falls back into the dark matter 
potential well of the main component. The infalling material 
expands adiabaticaUy as it fills up the available space at the 
center of the main cluster, resulting in much lower gas den- 
sity and temperature than the original ones. The subcluster 
core follows the shock and falls towards the center of the dark 
matter potential well of the subcluster. About 1 Gyr after the 
first core passage (t = 2 Gyr) the subcluster core is clearly sep- 
arated from the main cluster core, as found by Zu Hone et al. 
(2009b). 

In the present study, we shall focus on global physical 
properties of the cluster system, such as X-ray and Sunyaev 
Zel'dovich effect (SZE) observables, to be measured during 
the collision. We derive and compare values of the X-ray 
luminosity Lx and the integrated Comptonization parameter 
Y oc JydVnT, or the low-frequency integrated SZE flux, at 



times before and after the collision. 

We choose the direction of projection such that the centers 
of the two post-collision clusters are aligned along the projec- 
tion axis, to be consistent with the high-resolution X-ray ob- 
servations (Ota et al. 2004; Zhang et al. 2005). For the X-ray 
luminosity we assume Lx oc ex oc n^VT, with ex being the X- 
ray emissivity. The quantities Lx and Y are measured within a 
fixed aperture radius of rA^soo- We find that about t = 2{3) Gyr 
after the collision, the X-ray luminosity Lx and the integrated 
Comptonization parameter Y within rjoo are reduced to about 
11%(19%) and 33%(47%) of their respective initial values, 
reveahng substantial reductions in projected X-ray and SZE 
observables. 

In Figure 24 we show the projected X-ray spectroscopic- 
like temperature profiles (Mazzotta et al. 2004) for our simu- 
lated head-on cluster merger at three different times, namely 
f = 0,2, and 3 Gyr. Here the spectroscopic-like measure Tsi = 
JydVWsiT/ JydVWsi with the weight W^i = n^T-^l^ provides 
a good match to the temperature derived from X-ray spectro- 
scopic data (Mazzotta et al. 2004). We find that the projected 
X-ray temperatures Tx at later times, f = 2 - 3 Gyr after the 
beginning of the collision, are dramatically decreased by fac- 
tors of about 2, showing fairly flat radial profiles and thus re- 
sembling isothermal profiles as found by Chandra and XMM- 
Newton observations (Ota et al. 2004; Zhang et al. 2005). We 
note that the projected X-ray temperature in post-shock states 
is biased low due to the presence of low-temperature compo- 
nents in an adiabatic expansion process. The amplitudes of 
these post-shock Tx profiles are in good agreement with the 
observed average X-ray temperature Tx = 4.47^^ 54 keV (90% 
CL) from the Chandra observations of Ota et al. (2004), as 
compared to the values Tx = 2.6-2.8keV obtained by Zu 
Hone et al. (2009b). Since both the simulations presented 
here and in Zu Hone et al. (2009b) adopted the same mass 
ratio, impact parameter, and relative velocity as initial con- 
ditions, it is likely that the different post-shock temperatures 
at f = 2 - 3 Gyr are due to the different total masses of the 
two cluster components adopted. We plan to explore in more 
details the parameter space of our simulation, using different 
mass ratios, impact parameters, etc., in order to constrain the 
parameter space of possible solutions. 

8. SUMMARY 

In this paper, we have presented a joint weak and strong 
lensing analysis of the rich, but X-ray faint, cluster of galaxies 
C10024-H1654 at z = 0.395 based on wide-field Subaru BR^z' 
imaging combined with detailed strong lensing information 
obtained from deep HST/ACS/N1C3 observations (Zitrin et al. 
2009b). 

The deep Subaru three-band photometry, in conjunction 
with our weak lensing dilution techniques (Medezinski et al. 
2007; Umetsu & Broadhurst 2008; Medezinski et al. 2009), 
allows for a secure selection of distant blue and red back- 
ground populations free from contamination of unlensed 
galaxies, providing a greater lensing depth than achievable in 
the standard color-magnitude selection method. 

Our non-parametric mass reconstruction from a full lens- 
ing analysis of joint Subaru and ACS/N1C3 observations re- 
veals a continuously steepening density profile over a wide 
radial range from 40 to 2300 kpc with only a minor contri- 
bution, SM/M ~ 10% in the mass, from known substructure 
(Czoske et al. 2002; Kneib et al. 2003) at a projected distance 
of ~ 700 kpc /z"^ The cluster light profile closely resem- 
bles the mass profile, and our model-independent M/Lg pro- 
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file shows an overall flat behavior with a mean of {M/Lr) ~ 
2?>Qh{M / Lr)q, in contrast to centrally peaked M /L profiles 
found for other massive, relaxed clusters, and exhibits a mild 
declining trend with increasing radius at cluster outskirts, 
r ^ 0.6rvir. We found that the projected mass distribution 
for the entire cluster can be well fitted with a single NFW 
profile with virial mass, Mvir = (1.2 ±0.2) x IO'^Mq/z"', in 
agreement with that obtained from a model independent ap- 
proach (§6.2), and with those from recent lensing observa- 
tions (Hoekstra 2007; Zitrin et al. 2009b), but in apparent dis- 
agreement with X-ray hydrostatic mass estimates (Ota et al. 
2004; Zhang et al. 2005), and with caustic mass estimates at 
cluster outskirts (Diaferio et al. 2005). 

Careful examination and interpretation of X-ray and dy- 
namical data (Czoske et al. 2002; Ota et al. 2004; Jee et al. 
2007), based on detailed high-resolution cluster collision 
simulations (Czoske et al. 2002; Zu Hone et al. 2009a,b), 
strongly suggest that this cluster system is in a post colli- 
sion state, which we have shown is consistent with our well- 
defined mass profile for a major merger of two similar-mass 
clusters occurring along the line of sight, viewed approxi- 
mately 2-3 Gyr after impact (Zu Hone et al. 2009b) when 
the gravitational potential has had time to relax in the cen- 
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Table 1 . Properties of the Galaxy Cluster Cl0024+ 1 654 



Parameter Value 



Optical center position (J2000.0) 




R.A 


00:26:35.69 


Decl 


+17:09:43.12 


Redshift 


0.395 


X-ray temperature (keV) 


4.47!^ g (90% CL) 


Einstein radius (") 


30 (at z= 1.675) 



References. — [1] Czoske et al. (2002); [2] Ota et al. (2004); [3] Zitrin et al. (2009b). 

Note. — The optical cluster center is defined as the center of the central bright elliptical galaxy, or the galaxy 374 in the spectroscopic catalog of Reference 
[1]. Units of right ascension are hours, minutes, and seconds, and units of decUnation are degrees, arcminutes, and arcseconds. The average X-ray temperature 
Tx is taken from Reference [2]. The Einstein radius 6e for a background source atz = 1.675 is constrained by detailed strong lens modeling by Reference [3]. 



Table 2. Subaru Suprime-Cam data 



Filter 


Exposure time" 


Seeing'' 








(arcsec) 


(AB mag) 


B 


3 X 1200s 


1.27 


27.1 


Rc 


11 X 480s 


0.80 


26.8 


z' 


5x300s+3x60s 


0.82 


25.3 



"Total exposure time in units of s. 

''Seeing FWHM in the final co-added image. 

"Limiting magnitude for a 3it detection within a 2" aperture. 
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Table 3 . Weak lensing galaxy samples 



Sample name 


Magnitude limit" 


N 










{^2)g 




(AB mag) 




(arcmin^^) 










Red 


21.0 < z' < 25.5 


8676 


10.9 0.447 


1.14 ±0.09 


1.09 ±0.07 


0.56 ±0.02 


0.32 


Green 


17.4 < z' < 25.5 


1655 


2.1 0.396 


0.46 ±0.02 


0.45 ±0.02 


0.11 ±0.01 


0.03 


Blue 


23.0 < z' < 25.5 


5004 


6.3 0.463 


1.81 ±0.17 


1 75+<''*'' 


0.68 ±0.04 


0.47 


Blue+red 




13680 


17.2 0.453 


1.31 ±0.09 


• -0.30 
1 9040.16 
^•^^-0.14 


0.61 ±0.03 


0.38 



^Magnitude limits for the galaxy sample. 

''Mean surface number density of background galaxies in the C10024+1654 field. 
'^Mean rms error for the shear estimate per galaxy, Tfg = (ct^)'/^ 

''Mean redshift of the background sample derived with the COSMOS photometric catalog. 

"Effective source redshift corresponding to the mean depth of the COSMOS background sample. 

'Distance ratio averaged over the COSMOS redshift distribution of the background sample, (/3) = {D^jDs). 

^Distance ratio squared averaged over the COSMOS redshift distribution of the background sample, {0^) = {(D^s/Ds)^). 



Table 4. Summary of the best-fitting NFW/gNFW parameters 



Method 


WL" 


SL" E.R.'= (6>io,0up)'' 


a" 












(') 




(IO'SMq//,) 




(") 



Strong lensing 




IDk 




(0.17,0.79) 


1 


1 s;9+0,39 
i 27 


7 ^+1.4 


0.71/14 


34+13 

■^^-10 


Tangential shear 


lDg+ 






(0.55, 16) 


1 


1 14+0.22 
-QlS 




2.2/8 


37+11 
-IQ 








Yes 


(0.50, 16) 


1 


1 i9-(<).23 




4.9/9 


32+10 


Cc -Statistic 


IDk 
IDk 




Yes 


(0.55, 12) 
(0.50, 12) 


1 
1 


q 7+14.6 

8.6!|:' 


3.6/8 
3.6/9 


33+25 
■^^-20 
30+11 




IDk 


IDk 




(0.17,12) 


1 


^■^•'-0.15 


-1+1.4 
'•"'-1.2 


4.1/20 


33+f 




IDk 


IDk 




(0.17,12) 


04+0-93 




9.8 ±1.2 


3.8/20 


32i 



"Type of Subaru weak lensing (WL) data. 

''Fitting with or without the inner k profile derived from the strong lensing (SL) analysis of Zitrin et al. (2009b). 
"^Fitting with or without the inner Einstein-radius (E.R.) constraint, 9^ = 30" at z = 1.675. A 10% error is assumed for 9^. 
''Lower and upper radial limits of lensing constraints used for fitting. 

'^Central cusp slope and its 68.3% confidence interval including the uncertainty in the source redshift calibration. For NEW models, a is fixed at a = 1. 
'Virial mass and its 68.3% confidence interval including the uncertainty in the source redshift calibration. 

^Halo concentration, c-2 = Cvir/(2— a), and its 68.3% confidence interval including the uncertainty in the source redshift calibration. 
''Einstein radius in units of arcseconds for a background source at z = 1.675 as predicted by the best-fit NFW model. 



Table 5 . Best-fit parameters for the two-component NFW lens model 



Halo component 


Profile 


Mvir 


'^vir 


AR.A. 


ADec. 










(') 


(') 


Central 


NFW 


1.11±0.18 


8.1 ±1.2 


0.0 


0.0 


Northwest 


tNFW 


0.128±0.051 


5.0±3.5 


-2.3 


2.5 



Note. — Shown are the best-fit parameters and their Icr errors for the two-component NFW lens model (§5.2) derived from a joint fit to the inner k profile 
from strong lensing and the outer Subaru two-dimensional distortion data. The quoted errors include the systematic uncertainty in the redshift distribution of 
the background galaxies (Table 3). The central-halo centroid is fixed at the dark-matter center of mass (§4.2). The northwest-halo centroid is fixed at the peak 
position of the northwest galaxy clump in the surface number density distribution of BRcz'-selected cluster galaxies. We have applied a prior to the central halo 
to represent the constraints from inner strong-lensing information at 9 ^ 48" . A truncated form of the NFW profile (tNFW) is used for describing the projected 
lensing fields of the northwest halo. The resulting value is — 36761 with 23306 degrees of freedom. 



Table 6. Three-dimensional cluster mass from a deprojection analysis 



A Ma fA 

(lOi'Mo/A) (Mpc//i) 

"2500 0.46!<'<'3 

500 72+*' i' 93+0 07 
200 1.17!|j-| 

Virial 1.21!0:^| l-75i|:?5 



22 Mass Structure of C10024+1654 from a Full Lensing Analysis 




Fig. 1. — Subaru BRz' pseudo-color image of the central 8' X 8' region of the galaxy cluster C10024+1654 at z = 0.395. The side length of the field is 
1.8 Mpc /i"' at the cluster redshift. 
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Fig. 2. — The quadrupole PSF anisotropy field as measured from stellar ellipticities before and after the PSF anisotropy correction. The left panel shows the 
raw ellipticity field of stellar objects, and the right panel shows the residual ellipticity field after the PSF anisotropy correction. The orientation of the sticks 
indicates the position angle of the major axis of stellar ellipticity, whereas the length is proportional to the modulus of stellar ellipticity. A stick with the length 
of 5% ellipticity is indicated in the top right of the right panel. 
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Fig. 3. — Stellar elhpticity distributions before and after the PSF anisotropy correction. The left panel shows the raw ellipticity components (e* jfij) of stellar 
objects, and the right panel shows the residual eUipticity components {Se* , <5eJ ) after the PSF anisotropy correction. 
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Fig. 4. — Mean statistical weight Ug as measured from the magnitude-selected galaxy sample, shown as a function of the object Gaussian size rg (left) and of 
the Rc magnitude (right). In each of the panels, the statistical weight Ug is normalized to unity in the first bin. 
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Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 5. — Sample selection in the color-color diagram of C10024+1654, displaying the green sample (green points), comprising mostly cluster member galaxies, 
and the red (red points) and blue (blue points) samples, comprising of background galaxies. The galaxies that we identify as predominantly foreground lie in 
between the cluster and background galaxies are marked in magenta. 
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Fig. 6. — Left: Distribution of mean projected distances from the cluster center for all galaxies in the C10024+1654 field, displayed in color-color (CC) space. 
The bluer colors imply smaller mean radii, hence correspond to the location of the cluster in CC space. The black box marks the boundaries of the green galaxy 
sample we select which conservatively includes all cluster members. Right: Number density of galaxies in CC space. Several distinct density peaks are shown to 
be different galaxy populations - the reddest peak in the upper-right corner of the plot depicts the overdensity of cluster galaxies, whose colors are lying on the 
red sequence; the middle peak with colors bluer than the cluster shows the overdensity of foreground galaxies; the remaining peaks in the bottom part (bluest in 
B—Rc) can be demonstrated to comprise of blue and red (left and right, respectively) background galaxies. 



Umetsu et al. 



25 



CN - 



1 1 — I 1 — I 1 — I 1 — I 1 — I 1 1 — I 1 — I 1 — I r 



r 

COSMOS red 

COSMOS green 

COSMOS blue 



'.ri 



JiL. 



12 3 4 

COSMOS photometric redshift, z 



Fig. 7. — Photometric redshift distributions N(z) of BR^z' selected galaxy samples in the 2-deg^ COSMOS field. The green (green dotted line), red (red solid 
line), and blue (blue dotted-dashed line) galaxy samples are selected according to color-color/magnitude limits used in the weak-lensing analysis of C10024+1654. 
Each distribution is normalized by / dzN(z) - 1. 
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Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 8. — Gravitational reduced-sliear field in C10024+1654 obtained from shape distortions of the blue+red background galaxies, smoothed with a Gaussian 
with FWHM = 1.4' for visualization purposes. A stick with a length of 10% shear is indicated in the top right corner. The filled circle indicates the FWHM of 
the Gaussian. The coordinate origin is at the center of the cD galaxy. 




Fig. 9. — Comparison of the surface mass density field and the cluster galaxy distributions in C10024+1654. Left: Dimensionless surface mass density field, 
or the lensing convergence k(6) = Em(0)/Sciit, reconstructed from Subaru distortion data. Middle: Observed surface number density distribution S;,(0) of 
BScz'-selected green galaxies (17.3 < z' < 25.5 AB mag), representing unlensed cluster member galaxies. Right: Observed Rc-band surface luminosity density 
distribution S/(0) of the same cluster membership. The solid circle in each panel indicates the cluster virial radius of Syir — 8', or ry„ ~ 1 .8 Mpc /r' at the cluster 
redshift of z = 0.395. All images are smoothed with a circular' Gaussian of FWHM 1.4'. Also overlaid on the T,„(0) and S;(0) maps are the k(0) field shown in 
the left panel, given in units of 2cr reconstruction eiTor from the lowest contour level of 3(t. The field size is 26' X 22' . North is to the top, east to the left. 
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Fig. 10. — Azimuthally-averaged radial profiles of the tangential reduced shear g+ (upper panel) and the 45° rotated (x) component gx (lower panel) for 
our red (triangles), blue (circles), and green (crosses), galaxy samples. The error bars represent 68.3% confidence intervals estimated by bootstrap resampling 
techniques. The red and blue populations show a very similar form of the radial distortion profile which declines smoothly from the cluster center, remaining 
positive to the limit of our data, 6,i,a,K = 16'. For all of the samples, the X -component is consistent with a null signal at all radii, indicating the reliability of our 
distortion analysis. 




I I I I I d 

0.5 1 2 5 10 20 

9 [arcmin] 



Fig. 11. — Azimuthally-averaged radial profile of the tangential reduced shear g+ (upper panel) as measured from Subaru distortion data of our com- 
posite blue-H'ed background sample. The Einstein radius constraint (triangle) of 9^ = 30" ±3" (at Zs = 1.675), determined from multiply lensed images in 
HST/ACS/NIC3 observations (Zitrin et al. 2009b), is translated to the corresponding depth of the Subaru blue-Hred background sample (Table 3), using the best-fit 
NFW model to the Subaru and ACS/NIC3 data (see Figure 14), and added to the distortion profile (g+ = 1), marking the point of maximum distortion. The solid 
curve shows the best-fit NFW profile for the Subaru g+ measurements. The dashed curve shows the NFW profile from a joint fit to the inner Einstein-radius 
constraint and the outer Subaru g+ profile. Shown in the bottom panel is the 45°-rotated gx component. 



28 



Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 12. — Radial profile of lensing convergence k(9) = T,i„(9)/Y,^,\i as reconstructed from Subaru distortion data. The open squares show the results from 
a non-linear extension of aperture mass densitometry based on azimuthally-averaged tangential distortion measurements. The error bars are correlated. The 
gray-shaded region represents the Einstein-radius constraint on the mean interior convergence, k(< = 1, translated into the corresponding depth of the Subaru 
blue-Hred background sample. The dashed and soHd curves are the best-fitting NFW profiles from the Subaru k data with and without the inner Einstein radius 
constraint combined, respectively. 
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Fig. 13. — Number-count profile of BScz'-selected red galaxies {squares) in the background of C10024-I-1654. The triangles show the counts without the mask 
correction due to cluster members and bright foreground objects. The gray-filled region represents the 68.3% confidence bounds for the predicted count depletion 
curve from an NFW model constrained by our Subaru distortion analysis, demonstrating clear consistency between these two independent lensing observables. 
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Fig. 14. — Radial surface mass density profile of the galaxy cluster C10024+1654 over a wide range of radius from 40 to 2300 kpc /i~ reconstructed from our 
joint weak and strong lensing analysis of Subaru and ACS/N1C3 observations. All of the radial profiles are scaled to a fiducial redshift of = 1. The squares 
represent our Subaru results (this work) from a one-dimensional reconstruction by a non-Unear extension of aperture mass densitometry based on azimuthally- 
averaged tangential distortion measurements. The error bars are correlated. The triangles represent the inner k profile derived from the strong-lensing analysis 
of Zitrin et al. (2009b) based on 33 multiply-lensed images, spread fairly evenly over the central region, W ^jd ^ 48". The thick {black) and thin (gray) solid 
curves show the best-fitting NFW and gNFW profiles, respectively, from our full lensing analysis of ACS/N1C3 and Subaru observations. The best-fitting NFW 
profile from the inner k profile from the ACS/N1C3 observations is also shown as a {blue) dashed curve. The Subai'u weak-lensing constraint at the innermost 
radius is fully consistent with the strong-lensing information within the measurement uncertainty. 
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Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 15. — Joint constraints on the NFW model parameters (cvir, ^vir) for C10024+1654 derived from the radial profile of lensing convergence, = 
S,„(6)/S(.ri( (see Figure 14). The open square shows the best-fit set of the NFW model parameters, and the dashed contours show the 68.3%, 95.4%, and 
99.7% confidence levels in the cvir-Wvir plane. The solid contours show the same confidence levels, but for an NFW fit to the inner (10" 48") k profile 

constrained by deep ACS/N1C3 observations of Zitrin et al. (2009b). The corresponding best-fit set of (fyirj^vir) is shown by the open triangle. Also shown by 
the filled gray ai'eas ai'e the same confidence areas, but obtained for a joint fit to the combined ACS/N1C3 and Subaru «;(9) data. The corresponding best-fit set of 
(Cvir^^vir) is shown by the cross. 



LO 



0) 

o 



ft 

M 
O 

Id 

Pi 

CD 



LO 

d 



_L 




I 




0.5 1 1.5 2 



10 

c o 



15 



20 



Fig. 16. — Constraint on the the gNFW model pai'ameters, namely the central cusp slope a, the halo virial mass Myi,, and the halo concentration c_2 = 
Cvir/(2 — a), when all of them are allowed to vary, derived from the full lensing profile of CL0024-H1654 shown in Figure 14. The left and right panels show the 
two-dimensional marginalized constraints on {M^\^,a) and (c_2,a), respectively. In each panel of the figure, the contours show the 68.3%, 95.4%, and 99.7% 
confidence levels, and the cross indicates the best-fit model parameters. 
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Fig. 17. — Joint constraints on tlie NFW model parameters (Cvir,Mvir) derived from Subaru distortion data. The dashed contours show the 68.3%, 95.4%, 
and 99.7% confidence levels in the Cyii-Mvi, plane, estimated from Ax^ = "Xmin = 2.3, 6.17, and 11.8, respectively, for the tangential distortion (g+) profile 
(Figure 11). Also shown by the filled gray areas are the same confidence regions, but obtained for a joint fit to the outer Subaru g+ profile and the strong-lensing 
constraint on the location of the Einstein radius, = 30" ± 3" at Zs = 1.675. The open triangle shows the best-fit set of the NFW model parameters, (Cvir.Mvi, ), 
for the Subaru g+ results. The cross shows the best-fit set of (Cvji,Myii) for the combined Subaru g+ and Einstein-radius constraints. The solid curve shows the 
NFW c„ir-Mvir relation for 6»e = 30" at Zs = 1.675. 
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Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 18. — The i?c-band surface luminosity density profile S;(0) {squares) of B/?cz' -selected cluster member galaxies in C10024-I-1654. The gray-filled region 
represents the 68.3% confidence bounds for the surface mass density profile reconstructed from our joint weak and strong lensing analysis (Figure 14), 

converted into a luminosity density assuming a constant mass-to-light ratio of (M/L^)(< 3') ~ 230/i(M/Lr)q. Also shown with filled circles is the same cluster 
luminosity density profile, but corrected for the presence of the northwest clump located at a projected distance of 9 ~ 3'. 
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Fig. 19. — Model-independent radial profile of the differential mass-to-light ratio &M{Q) j &Lr(Q) = S,„(9)/E;(6) derived using the surface mass density profile 
Sm(9) (see Figure 12) from the Subaru weak-lensing analysis and the X'-coirected Sc-band surface luminosity density profile S;(6) measured from BRcz'-selected 
cluster member galaxies. The error bars represent 68.3% confidence limits. 
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Fig. 20. — Comparison of tangential distortion (g+) profiles from different lensing studies. No correction is applied to the mean depth of lensing observations 
between different studies. The hatched region shows the 68.3% confidence interval for the ACS strong-lensing results based on the azimuthally-averaged k profile 
of Zitrin et al. (2009b), scaled to a source redshift of Zs = 1.3, roughly matching the mean depth of the Subaru blue+red background sample (see Table 3). Also 
shown with the solid curve is the best-fit NFW model (zs = 13) from the ACS strong lensing constraints. The squares show our Subaru g+ profile as measured 
from BRcz'-selected blue+red background galaxies. A simple extrapolation of the ACS/NIC3-derived NFW profile (solid curx'e) fits well with the outer Subaru 
distortion information {squares) over a wide range of radius, but somewhat overpredicts the distortion profile at 3' 9 ,^ 5'. The triangles show the g+ profile 
from the weak lensing analysis of Jee et al. (2007, Figure 13) based on deep 6-passband HST/ACS images, where the data are limited to the positive parity region 
(9 ^ 40") in this comparison. The open circles show the g+ profile at 0.9' ^0 ^ 11' from the weak lensing analysis of Kneib et al. (2003, Figure 7) based on a 
sparse-sampled mosaic of 2-band (F450W, F814W) WFPC2 observations. The flattened slope of Kneib et al. g+ profile (zs = 1.15 it 0.3) at ^ 2' is likely due 
to contamination of the weak lensing signal by unlensed cluster member galaxies. Note that Kneib et al. (2003) made a correction for this dilution effect in their 
two-dimensional lens modelling. 
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Fig. 2 1 . — Comparison of projected mass density (k) profiles of C10024+1654 from different lensing studies. All of the radial profiles are scaled to a fiducial 
redshift of = 1. The squares represent our model-independent k(8) profile (this work) reconstructed from Subai'u distortion data (Figure 12). The error bars 
are correlated. The triangles represent the results from the strong-lensing analysis of Zitrin et al. (2009b) based on deep ACS/NIC3 images. The solid curve 
shows the best-fitting NFW profile from our full lensing analysis of ACS/NIC3 and Subaru observations covering a wide range of radii, R ~ [40, 2300] kpc /i"' . 
The dashed curve shows the k profile (0 ^ 1.8') reconstructed from the ACS weak and strong lensing analysis of Jee et al. (2007). The dotted curve shows 
the two-component lens mass model constrained from the WFPC2 weak lensing measurements of Kneib et al. (2003) combined with the inner strong-lensing 
constraint on the Einstein radius. 
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Mass Structure of C10024+1654 from a Full Lensing Analysis 
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Fig. 22. — Comparison of cumulative projected mass profiles M2d(< Q) of the galaxy cluster C10024+1654 from different lensing studies. The solid curve 
(this work) represents the results from our joint weak and strong lensing analysis of deep Subaru and ACS/NIC3 observations over a wide radial range from 
40kpc/i"' to 2300kpc/i"'. The gray-shaded area shows the 68.3% confidence interval at each radius estimated from a Monte-Carlo error analysis taking into 
account the error covariance matrix of our full lensing constraints. Shown with a dotted vertical line is the maximum radius of multiply-lensed images used in the 
strong lensing analysis of Zitrin et al. (2009b). The triangles represent the results from the ACS weak lensing analysis of Jee et al. (2007, Figure 12). The dashed 
curve shows the best-fit two-component lens model of Kneib et al. (2003), constrained from their WFPC2 weak lensing measurements combined with the inner 
Einstein radius constraint, taking into account the contribution of both the central and northwest clumps identified in projection space. The three mass profiles 
from different lensing work are in good agreement at the Einstein radius 6^ ~ 30" (as indicated by the arrow) for the 5-image system (z, = 1.675), by which each 



profile is normalized. The projected mass profile of Jee et al. (2007) is in an excellent agreement with our joint mass profile out to ( 
projection space), but increasingly exceeds our profile 3X6 ^ 1.1' out to the liinit of their data. 
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Fig. 23. — Comparison of cuinulative spherical inass profiles M3d(< 9) of the galaxy cluster C10024-I-1654 from different studies. The solid curve (this work) 
represents the deprojected spherical mass profile from our full lensing analysis of Subaru and ACS/NIC3 data assuming spherical symmetry for the cluster system. 
The gray shaded area shows the 68.3% confidence interval at each radius estimated from a Monte-Carlo error analysis taking into account the error covariance 
inatrix of our joint Subaru and ACS/NIC3 lensing constraints. Shown with a dashed curve is the spherical NFW model of Kneib et al. (2003) for the central 
component alone (see also Figure 22). The triangles represent the results from the dynamical analysis of Diaferio et al. (2005) based on the caustic method. The 
dotted curve shows the best-fitting NFW profile from the Chandra X-ray analysis of Ota et al. (2004) obtained assuming the hydrostatic equilibrium. Note that 
the maximum limit of the Chandra data of Ota et al. (2004) is about 4' (900 kpc /?"'). Our joint mass profile is only marginally consistent with the dynamical 
results of Diaferio et al. (2005) out to (? ~ 3.5', but increasingly exceeds the caustic-based profile at large radii. A significant increase in M3|3(< r) is seen in our 
joint mass profile at 3.5' ^9^5', indicating an additional, extended mass component in the outer radius. 
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Fig. 24. — Projected spectroscopic-like temperature profiles created from mock X-ray observations of our simulated head-on 2:1 mass-ratio cluster merger, 
shown at three different times, ? = 0, 2, and 3 Gyr after the beginning of the collision (open squares, crosses, and open triangles, respectively). The vertical en'or 
bars represent the dispersion due to azimuthal averaging. The hatched region shows the 90% confidence interval for the average X-ray temperature, Tx = 4.47;*^ 
from Chandra X-ray observations of Ota et al. (2004). 



